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COMPUTER  PROGRAM  FOR  MICROWAVE 
GaAs  MESFET  MODELING 


INTRODUCTION 

There  has  long  been  a  need  for  a  simple,  inexpensive,  fairly  accurate,  easy-to-run  computer  pro¬ 
gram  for  calculating  the  equivalent  circuit  elements  of  MESFETs  (metal-semiconductor  field  effect 
transistors).  Microwave  circuit  designers  use  an  equivalent  circuit  representation  for  the  transistor, 
similar  to  that  shown  in  Fig.  1,  which  will  be  discussed  in  detail  later.  They  also  sometimes  use  S- 
parameters  to  represent  the  active  device,  but  they  can  be  calculated  directly  from  the  equivalent  cir¬ 
cuit.  A  problem  with  all  the  analytical  MESFET  models  proposed  to  date  is  that  they  do  not  give 
answers  which  agree  well  with  measured  data  over  all  bias  conditions.  Provisions  have  been  made,  in 
the  computer  program  FETREN,  presented  in  this  report,  to  improve  on  agreement  with  measured 
results  over  the  full  range  of  bias  conditions  (or  over  the  full  range  of  large  signal  ’’swing"  conditions). 


Fig  I  —  Equivalent  circuit  represeniation  of  a  GaAs  MESFET.  The  most  significant 
bias-dependent  device  elements  are  C  C'/*,  ,em.  K/*.  and  Rq 

Certain  new  terms  are  added  to  the  model  employed  in  this  program  to  permit  adjustment  of  com¬ 
puted  results  to  measured  values— mostly  in  a  self-consistent  and  plausible  way,  rather  than  by  the  use 
of  arbitrary  "fitting"  constants.  In  areas  where  the  analytical  theory  is  not  exact  enough,  or  where  the 
knowledge  of  the  field  behavior  is  not  well  enough  known,  computational  methods  of  interpolation  are 
employed  to  obtain  estimates  of  the  circuit  element  values  between  known  values.  For  example,  this 
approach  is  required  to  get  reasonable  answers  when  the  drain  voltage  is  somewhere  between  low  values 
for  which  no  Gunn  domain  can  exist,  and  high  values  for  which  the  domain  is  "fully"  formed. 

Other  such  departures  from  earlier,  less  practical,  models  will  be  described  subsequently.  The 
philosophy  applied  in  the  development  of  the  model  used  for  the  computer  program  FETREN  is:  (1> 
apply  whatever  analytical  theory  gives  the  best  agreement  with  measured  results,  (2)  do  not  be  disap¬ 
pointed  with  less  than  perfect  agreement,  and  (.1)  keep  the  program  simple  and  inexpensive  to  use 
The  program  listing  is  included  as  Appendix  A,  and  the  flow  diagram  is  included  as  Appendix  B. 
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Table  1  describes  the  symbolic  notations  used  in  the  text  in  the  order  in  which  they  appear.  A 
notation  list  for  the  computer  program  is  included  with  the  program  listing  in  Appendix  A. 

Table  I  —  Text  Notation 
8  Conducting  channel  thickness  in  Region  II 

Potential 

p  Volume  charge  density 

« ,  Relative  dielectric  constant 

€  o  Permittivity  of  free  space 

«  Product  of  «0  and 

Vos,  Longitudinal  voltage  drop  across  Region  1 

•'av,  Longitudinal  voltage  drop  across  Region  II 

Longitudinal  voltage  drop  along  the  conducting  channel 
from  the  drain  end  of  the  gate  to  the  source  end  of  the  gate 

y  Vertical  distance  measured  from  bottom  of  active  layer 

b  Vertical  thickness  of  conducting  channel 

a  Total  active  layer  thickness 

Vpo  Pinch-off  voltage 

q  Electronic  charge 

N0  Volume  donor  impurity  concentration 

Pn  Low  field  mobility 

IT'  Gale  width 

x  Horizontal  distance  measured  from  source  edge  of  gate 

fi)s  Drain-to-source  current 

L  Gate  length 

L  |  Length  of  Region  1  under  gate 

L2  Length  of  Region  11  under  gale 

/,  A  specific  analytical  function 

V'(l s  Gate-to-source  voltage 

Vm  Metal-to-semiconductor  gate  schottky  barrier  built  in 

voltage,  taken  as  a  positive  voltage 

b |  Conducting  channel  thickness  at  first  point  right  of 

source  where  the  field  has  increased  to  Es 

L's  Longitudinal  field  strength  which  just  causes  carrier 

velocity  saturation 

p  Square  root  of  the  ratio  of:  (1)  the  voltage  at  the  junction  point 

between  Regions  I  and  II,  relative  to  the  gale  electrode, 
to  (2)  the  pinch-off  voltage 
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Table  1  —  Text  Notation  (Continued) 

s  Square  root  of  the  ratio  of:  (I)  the  voltage  at  the  source  end  of 

the  gate,  relative  to  the  gate  electrode,  to  (2)  the  pinch-off  voltage 

£,„„v  Longitudinal  electric  field 

R s  Sum  of  bulk  semiconductor  resistance  and  contact  resistance 

from  source  edge  of  gate  to  the  source  metal 

Ri,  Sum  of  bulk  semiconductor  resistance  and  contact  resistance 
from  drain  edge  of  gate  to  the  drain  metal 

k’o.v  Voltage  drop  from  gate  metal  to  the  conducting  channel 
at  the  source  end  of  the  gate 
Voltage  drop  from  gate  metal  to  source  metal 
I'as  Voltage  drop  from  drain  metal  to  source  metal 

AK  Substrate  current  multiplying  constant 

Isib  A  current  said  to  be  passing  through  the  substrate  material, 

under  the  active  layer,  in  parallel  with  the  active  current  generator 

DR  1  Drain  shunt  resistance 

vM,  The  assumed  electron  saturated  velocity 

RTH  Device  overall  thermal  resistance 

gm  1  ransconductance 

fg  A  specific  analytical  function 

d  Square  root  of  the  ratio  of:  (1)  the  voltage  at  the  drain 

end  of  the  gate,  relative  to  the  gate  electrode,  when  there 
is  no  Region  II,  to  (2)  the  pinch-off  voltage 

r0  Intrinsic  transistor  output  resistance 

gwb  An  output  conductance  term  related  to  substrate  current 

rsub  The  reciprocal  of  g,ub 

R//v  The  charging  resistance  for  Qv 

Clbi  The  gate-to-source  capacitance 

Q  Total  charge 

Qg  Charge  in  depleted  area  under  gate 

Q,,  Charge  of  extended  depletion  area  toward  drain 

Qdom  Charge  of  domain 

Cx,0  Gate-to-source  capacitance  for  drain  voltage  near  zero 

C'g Gate-to-source  capacitance  Tor  large  drain  voltage 

nit  The  "characteristic"  doping  density,  which  depends  on 

the  shape  of  the  velocity  vs  electric  field  curve,  see  Ref.  12 

A  (A.  )  Complete  elliptic  integral  of  the  first  kind  in  terms  of 
the  modulus,  k 
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DL 

DGL 

HG 

SL 

SGL 

SDL 


Qsw 

Csw 

Qx 

RDOM 

CDOM 

n 

J 

t 

v 

D 

H'.h 

E,h 


Table  1  —  Text  Notation  (Continued) 

Drain  metal  length 

Distance  between  drain  and  gate  metal 
Gale  length,  same  as  L 
Source  metal  length 

Distance  between  source  and  gate  metal 

Distance  between  source  and  drain  metal  (SDL  =  SGL  +  WG  +  DGL) 
Radius  of  sidewall  capacitance  area 

Total  charge  depleted  from  source-end,  sidewall,  quarter  circle  depletion  area 

Sidewall  capacitance 

Same  as  Csv 

Domain  resistance 

Domain  capacitance 

Volume  electron  density 

Current  density 

Time 

Drift  velocity 

Diffusion  coefficient 

Threshold  depletion  layer  thickness 

Threshold  field  for  bulk  negative  differential  mobility 

Transit  time 


CONCEPTUAL  REVIEW 

Much  has  been  written  about  the  general  evolution  of  microwave  MESFETs  and  their  modeling. 
For  an  excellent  review  see  Liechti  [ll.  A  repeat  of  that  is  not  undertaken  here.  Rather,  a  review  of 
the  specific  concepts  used  in  the  development  of  the  computer  program  FETREN  and  its  associated 
model  is  given.  The  origins  of  all  the  related  concepts  are  referenced  and  discussed  if  they  were 
derived  from  earlier  literature,  and  are  explained  if  they  are  new  ideas  developed  by  the  authors  of  this 
report.  It  is  assumed  that  the  reader  is  familiar  with  the  subject  or  will  read  the  referenced  material  as 
required. 

The  majority  of  the  calculations  are  based  on  the  theory  of  Grebene  and  Ghandi  [2].  They 
expanded  on  the  work  of  Shockley  |3],  as  Shockley  himselT  suggested  it  might  be  done,  to  allow  two 
distinct  regions  under  the  gale.  Region  I,  nearer  the  source,  is  taken  to  have  constant  carrier  mobility; 
and  Region  11,  nearer  the  drain,  has  saturated  carrier  velocity.  Figure  2  shows  the  configuration,  with 
the  point  which  separates  Region  1  from  Region  11  being  free  to  move  longitudinally,  so  that  it  can 
position  itself  at  the  first  instance  of  longitudinal  electric  field  large  enough  to  produce  velocity  satura¬ 
tion.  The  longitudinal  electric  field  in  Region  I  is  always  least  at  the  source  end  and  is  larger  at  posi¬ 
tions  closer  to  the  drain  end.  The  vertical  conducting  channel  thickness  in  Region  11,  ft.  is  assumed 


4 


NRL  REPORT  8561 


SOURCE  DRAIN 


t  ig  2  —  Assumed  MESl'ET  configuration 
for  theoretical  analysis 


constant  over  the  whole  of  that  region.  The  method  of  solving  the  applicable  equations  for  Regions  I 
and  II,  separately,  and  applying  suitable  continuity  criteria  between  the  two  regions  is  given  in  Grebene 
and  Ghandi  [21  and  is  extremely  well  detailed  in  Pucel  et  al.  [4). 

From  a  conceptual  standpoint,  the  theory  or  model  which  is  deduced  from  Refs.  2,  3,  and  4  has 
some  problems,  even  though  it  may  never  be  greatly  improved  upon  for  small  signal  calculations  at 
moderate  drain  voltages  if  no  Gunn  domain  forms.  Some  obvious  problems  are:  (1)  no  computational 
provision  for  the  presence  of  a  Gunn  domain,  (2)  no  provision  for  adjustment  of  parameters  due  to 
heating  from  the  bias  power,  (3)  no  provision  for  the  effects  of  extension  of  the  depletion  region  longi¬ 
tudinally  beyond  the  gate,  (4)  no  specific  provision  for  a  gate  to  source  shunting  resistance  arising  from 
either  simple  "leakage"  or  from  carrier  fringing  around  the  extremely  high  field  region  when  a  Gunn 
domain  is  present,  (5)  no  provision  for  the  effect  of  the  incorrect  approximation  that  the  channel  thick¬ 
ness  in  Region  II  is  nearly  zero  at  all  times  (this  approximation  is  suggested  in  Refs.  3  and  4),  (6)  no 
way  of  knowing  the  effects  on  the  Region  II  calculations  for  a  drain  terminal  whose  physical  shape  is 
not  like  the  one  Shockley  assumed  in  his  original  formulation  (a  semicircle  whose  shape  is  approxi¬ 
mately  that  of  the  assumed  drain  equipotential  shape). 

It  is  fairly  safe  to  say  that  there  are  too  many  variables  to  expect  to  ever  get  near-perfect  agree¬ 
ment  between  any  model  and  measured  data.  If  the  agreement  needed  must  be  extremely  good,  it  is 
probable  that  some  fitting  to  measured  data  is  required. 

All  of  the  problems  mentioned  above  were  addressed  in  the  development  of  the  program 
FETREN.  The  way  they  were  handled  and  the  effects  of  considering  them  and  making  some  provisions 
for  them  can  be  seen  in  the  following  section  as  the  circuit  element  calculation  methods  are  described. 
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CIRCUIT  ELEMENT  CALCULATIONS 
The  Operating  Point 

The  operating  point  is  set  by  the  longitudinal  position  of  the  dividing  line  between  Region  1  and 
Region  11,  sueh  that  the  sum  of  the  calculated  voltage  drops  across  the  two  regions  equals  the  total 
applied  voltage.  In  Fig.  3,  l  /)S|  +  l  /)Si  =  I  /)s  .  From  the  outset,  FFTREN  has  been  set  up  to  work 
from  the  external  values  of  the  two  bias  voltages— that  is,  those  values  which  would  be  measurable  in 
the  laboratory.  Therefore,  in  the  initial  calculations  to  determine  the  operating  point,  the  "parasitic" 
resistances  must  be  included  They  are  calculated  between  lines  1359  and  1368  in  the  program,  the  list¬ 
ing  for  which  is  given  in  Appendix  A,  with  the  drain  and  source  contact  resistances  included  in  Rs  and 
Rn  The  calculation  method  is  taken  from  Fukui  (5], 


Fig.  3  -  MESFET  configuration  showing  internal 
and  external  voltages 


An  applicable  relationship  in  the  depletion  region  of  Region  1,  for  a  volume  charge  density  p,  is 
Poisson's  equation: 


V2<t>  =  - 


(1) 


where  <l>  is  poicntial,  <,  is  the  relative  dielectric  constant  of  the  intrinsic  material,  and  e0  is  the  permit¬ 
tivity  of  free  space  Then,  by  the  gradual  channel  approximation  [3],  the  potential  variation  in  the  y- 
direction  is  much  greater  than  that  longitudinally;  the  double  integral  of  Eq.  (1)  in  the  .v-dimension 
only,  with  appropriate  boundary  conditions,  yields  potentials  for  Region  I: 


<ji  = 


L>  y2  +  1! 
«  2  6 


hy  4 


ii  £i  _  fL 

f  2  e 


ab\ 


(2) 


where  n  is  the  total  active  layer  thickness  and  h  is  the  channel  thickness  as  a  function  of  horizontal 
position  The  potential  along  the  boundary  between  the  depletion  region  and  the  conducting  channel  in 
Region  I  is  given  from  Eq.  (2)  for  v  =  h. 
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<1>, 


2f 


2 

.  /> 

1  ^ 

1  — 

=  (/'O 

I - 

a 

a 

(3) 


where  t>  -  qbin,  with  q  =  electronic  charge  =  1.602  x  10  19  C;  Nn  =  volume  doping  density;  e  =  eti<,, 
with  t„  =  8.854  x  10  14  F/cm  and  e,  =  12.9  according  to  Ref.  6;  and  Fw  =  (qN„a2)/( 2c)  is  defined  as 
the  pinch-off  voltage. 

The  current  in  the  conducting  channel  of  Region  I  is  given  by  Ohm's  law.  such  that 

/tty,  =  <IP-i)Ni)bW  ~Hx'  ^ 

where  d<l>/t).v  is  the  longitudinal  electric  field,  W  is  the  total  gate  width  of  the  device,  and  nu  is  the  car¬ 
rier  low  field  mobility.  After  an  integration  over  the  length  Lx, 

q^NoWa  (J) 


h,s- 


I-  i 


where 


with 


and 


A  =  P2  ~  -  y  ( P 3  -  s}). 


p'=  <  -  f'o.v  +  VosjnVfo) 

^  as )/ 1  Vpo  )  • 


^  = 


(f 


(6) 


(7) 


(8) 


in  terms  of  intrinsic  voltages,  with  F#,  =  metal-lo-semiconductor  gate  Schottky  barrier  built-in  voltage, 
taken  as  a  positive  voltage.  Also,  Fw  is  taken  as  a  positive  voltage. 


In  Region  II,  just  at  the  boundary  with  Region  1,  the  current  may  be  calculated  as 

lnsn  =  WES, 


(9) 


where  £v  is  the  longitudinal  field  which  just  causes  carrier  velocity  saluration,  and  b ,  is  the  channel 
thickness  at  that  point.  Since,  as  shown  in  Ref.  4, 

b[  =  a  ( I  —  p),  (10) 

l/>s„  =  (•“/>)  (11) 

If  we  equate  F.qs.  (11)  and  (5),  to  require  current  continuity  at  the  junctions  between  Region  1  and  II, 
then 

VpoVO 


Es ( 1  -  p) 


(12) 


Then,  to  gel  voltage  drops  involved,  the  integral  of  the  longitudinal  electric  field  over  the  length 
of  Region  I  is  used  to  get  Vm . : 


/v-t-i  rx~L  i 

,-n  =  ^/w,  =  J(_0  - 


r  w  i 

</4>, 

J.v-0 

cix 

cix. 


(13) 
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Also,  in  Region  11,  the  potential  distribution  can  be  found  vs  ith  the  solution  of  Laplace's  equation 
in  two  dimensions: 


V-’  <l>n  =  0. 


( 141 


However,  an  exact  solution  of  this  for  the  general  case  is  difficult,  and  in  Refs.  2.  3,  and  4  the  same 
approximate  solution  is  used,  originally  proposed  by  Shockley  (2? i .  The  adequacy  of  the  approximate 
solution  requires  that  the  physical  shape  of  the  drain  conductor  be  "somewhat  like"  the  shape  of  the 
equipotential  lines  far  away  from  the  drain  This  is  not  the  case  in  planar  MESFETs.  so  in  MIR  I  \ 
some  liberty  is  taken  with  a  constant  related  to  this  portion  of  the  overall  solution.  The  constant  (T)R2 
in  the  program  FETREN,  at  line  1415,  may  be  varied  to  affect  the  amount  of  voltage  drop  assigned  to 
Region  II  (but  the  form  of  the  analytical  expression  is  maintained)  .Specifically. 


‘Wii  , 

- - ax, 

d.x 


(15) 


for  which  an  approximate  solution,  at  y  =  0,  is: 

<„.*,=  —  l-s  sinh 

7 T 


nU.  -  /.,) 
hi 


Also,  from  F.q.  (13)  in  Region  1,  as  shown  in  Pueel  el  al.  [4], 

I/m.  =  l/’f)l/r  ~  v'* 


(16) 


(17) 


So,  the  total  (internal)  voltage  drop  across  the  intrinsic  transistor,  alter  substituting  F.q.  (12)  into  Eq 
(16),  is: 


\ns  =  l,,()(/r  -  s’)  +  —  /:s  sinh 

7T 


77 

2  a 


m 


(1  -  p)Es 


(18) 


Figure  4,  a  simplified  transistor  schematic,  shows  that  the  accessible  (measurable)  drain-to-source 
voltage  is  the  sum  of  the  intrinsic  part,  F/)s  .  and  the  drops  across  Rs  and  Rn.  Note  also  that  the 
accessible  gate-to-source  voltage  is  sum  of  the  intrinsic  part.  Fav  and  the  drop  across  Rs.  Therefore 
the  accessible  drain-to-source  voltage  is 

I'm,,.  =  1  ,>s ...  +  /,. s<«v  +  «;,)■  <191 

Also,  the  accessible  gate-to-source  voltage  is 


E<.s 


=  1, 


+  //,s(  R  s)  ■ 


(20) 


In  the  program  FETREN,  if  there  is  a  Region  II,  the  determination  of  the  overall  operating  point  is 
determined  by  forcing,  through  iteration,  the  applied  F/)s  >|  to  equal  the  terms  ol  Eq  (10)  at  line 
number  520(1,  while  forcing  simultaneously  the  applied  •  liV  to  equal  the  terms  ol  Eq.  (20)  at  line 


6010  The  values  of  s,  />,  and  //)s  are  the  final  values  which  the  iterative  process  assigns  them.  The 
value  of  s  is  calculated  at  line  3100,  />  by  iterative  assignment  it  4700,  and  //)S,  indirectly,  at  5200  and 
repeated  explicitly  at  620(1.  Also,  within  this  iterative  loop  which  extends  approximately  from  line 
2200,  where  the  external  gate-to-source  voltage  is  set,  to  line  6040  where  the  viteration  loop  ends,  the 
value  of  /.  |  is  determined  at  line  4000  Because  the  solution  method  is  iterative,  l.\  may  sometimes 
overshoot  and  be  temporarily  assigned  a  value  larger  than  the  total  gale  length  /,.  even  though  the  final 
value  may  he  less  than  /.  Therefore  the  program,  at  line  4000.  allows  /.,  to  be  greater  than  I.  Alter 
final  convergence,  however,  if  the  final  value  for  l,\  is  larger  than  /..  as  tested  at  line  6100,  the  indica¬ 
tion  is  that  the  bias  levels  were  such  that  there  is  no  Region  If  For  this  case  it  is  then  necessary  to 
completely  re  solve  the  problem  by  using  the  equations  that  are  applicable  for  Region  I  only  1  he 
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tig  4  —  Simplified  transistor  schematic 


authors  derived  these  equations  independently  and  found  them  to  agree  with  those  developed  in  Refs  ' 
and  8.  those  references.  Ref.  4.  and  the  TTTRTN  listing  may  be  consulted  for  development  formula¬ 
tion.  I  bis  portion  of  I  TTRT.N  starts  at  line  6620.  and  the  operating  point  is  determined  hv  iterating 
the  current  and  finding  a  simultaneous  combination  of  s.  J,  and  //>.,.  for  which  the  applied  current 
equals  the  analytical  value: 


//)S  —  U</(4  l>',/>  i  *  HI 


v-  -  ~  U'  -  s') 


(2D 


where  cl  is  now  related  to  the  voltage  at  the  drain  edge  of  the  gate.  In  terms  of  intrinsic  bias  voltages. 

</•  -  <  *  /<,  -  <„s  +  r„s>./<  I  ,,,).  (22) 

analogous  to  Kq  (7).  In  summary,  then,  the  operating  point  is  prescribed  as  a  function  of  the  two  bias 
voltages,  as  outlined  above,  after  which  the  equivalent  circuit  parameters  and  the  actual  current  can  be 
computed. 


The  1-V  Curves 


The  current  calculated  in  the  preceding  section  is  an  initial  estimate  of  the  actual  current.  Other 
factors  must  be  considered  in  arriving  at  a  better  estimate  of  a  value  to  agree  with  measured  data.  If 
there  is  a  stationary  Gunn  domain  present,  a  substrate  current  resulting  from  electron  injection  into  the 
adjacent,  otherwise  assumed  semi-insulating,  substrate  may  result  |d|  The  arguments  in  Ref.  9  are  not 
very  strong  for  qualitative  accuracy.  Therefore,  in  TTTRTN.  a  constant  1^.  entered  at  line  1620  and 
used  at  lines  6206  or  6845,  is  included  to  make  the  substrate  current  eflect  variable  — but  the  lorm  of 
the  analytical  expression  is  maintained: 

/s,*«  (l,,s..„)  '  H.  <2J> 

This  shows  that  the  substrate  current,  in  the  presence  of  a  stationary  Gunn  domain,  is  proportional  to 
the  gate  width,  to  the  square  root  of  the  drain-to-source  voltage,  and  to  the  fourth  root  of  the  doping 
level  Tor  further  detail,  see  T.astman  and  Shur  [Of  (Substrate  injection  current  can  be  eliminated  by 
setting  I h  to  zero  at  line  1620.) 
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One  additional  term  is  added  to  the  currents  already  calculated  to  account  for  two  physical 
phenomena  not  otherwise  considered.  One  of  d'-se  phenomena  is  simple  resistive  leakage  across  the 
top  surface  of  the  transistor  or  through  ".  juustraie  whose  overall  conductivity  is  higher  than  desirable 
The  second  is  the  resistive  shunting  effect  of  the  doping  "tail,"  which  is  typically  disregarded.  The  dop¬ 
ing  tail,  in  Fig.  5,  is  the  crosshatched  area.  It  arises  because  of  the  use  01:  (I)  the  uniform  doping 
assumption  (adequate  for  almost  all  purposes),  creating  the  horizontal  dashed  line  at  the  doping  level  to 
be  used  in  the  analytical  expressions;  and  (2)  the  specific  active  layer  thickness  assumption,  creating  the 
vertical  dashed  line  at  the  thickness  to  be  used  in  the  calculations.  The  carriers  "left  over"  constitute  a 
shunting  resistor  in  parallel  with  whatever  other  resistive  shunting  terms  are  present.  In  FFTRFN,  pro¬ 
vision  is  made  for  a  fixed  resistance  to  shunt  the  intrinsic  device,  in  addition  to  the  possible  substrate 
current  around  a  Gunn  domain.  The  value  of  this  resistance,  assumed  proportional  to  gate  width,  is 
assigned  to  DR  I  at  lines  1420  and  1430,  and  its  contribution  to  the  drain  current  is  calculated  at  lines 
6216  (when  there  is  a  velocity  saturated  Region  ID  and  6862  (when  there  is  no  Region  ID. 


I  ig  '  —  I- sample  doping  profile  showing  doping  Tail” 


In  Rets  2  and  4  it  is  hinted  that  typically  the  finite  thickness  of  the  conducting  channel  in  Region 
II  may  be  disregarded  in  the  solution  of  Fq  (15)  However,  the  authors  carried  out  the  derivation  giv¬ 
ing  the  channel  the  actual  thickness  it  needs  to  carry  the  current  which  is  Rowing: 


where  v  is  the  assumed  electron  saturated  velocity  and  6  is  the  channel  thickness.  In  the  program 
this  appears  at  line  6225  Ihen  at  line  6235,  the  program  term  (  DR  1  is  calculated;  it  represents  the 
traction  ol  the  active  layer  thickness,  a.  which  is  to  be  used  in  the  next  iteration.  Notice  that  CDR  1 
was  preset  to  unity  at  line  1410,  outside  the  loop  which  determines  the  operating  point  and  the  drain 
current. 
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One  Iasi  ailllisiment  is  made  allot  the  lirsl  time  iht.uign  m  fe  i  mt  pmo:  dim:  vuiten:  limn 
H.ised  on  tile  exlcrn.ilh  applied  dram  voltage  an  I  l!)e  pielnttm  :  •  *  u,  i\,l  dim,  :he  power 

dissipation  is  ealeulaled  as  simph  the  product  ol  those  ,w..  tetnis  I  •  -  s  1  ■  •  •  i ..  n  poigi.nit  at  line 

<>’.'11  and  OXOX  I  hen  the  low  field  mohilm  value  is  adinsK  .1  ■■  n  pm  cm  i,se  icxuhmg  Iroin 

the  bias  power  dissipation,  as  unhealed  in  Appendix  <  at  lutes  t,.'4o  md  ox".'  Note  : i •  a t  a  tvpixul 

value  lor  the  thermal  testsianee.  Rill,  ol  a  \1I  SI  I  I  is  uKluded  m  the  pmei  mi  line  l-|so  x  sue 
jested  in  Appendix  and  that  inverse  proportion. dm  to  gate  width  is  assumed  t  I  herm.il 
modification  ol  mobility  may  he  eliminated  hv  setting  thermal  resistance  lo  /eto  at  I4'U  i  \  single 
lirsl  order"  correction  is  all  that  is  needed  for  the  Region  II  channel  thickness  and  lot  the  low-held 
mobility,  file  mobility  correction  for  bias  power  dissipation  is  primarily  responsible  for  the  downturn 
of  the  IV  curves  in  fig  C'2  of  Appendix  while  without  n  the  substrate  iniection  current,  around  the 
Ciunn  domain,  would  cause  the  curves  to  rise  w,th  increasing  drain  voltage 

Transconductance 


When  there  is  velocity  saturation,  or  when  there  is  a  Region  II  under  the  gate,  the  transconduc¬ 
tance  x,„  is  calculated  according  to  Slat/  et  al.  (10)  as: 


9  »  ,,s 


l  />s«  const 


otjH  „A/<  li  /:\ 

t  m 


(25) 


where 


[0  -  s)  cosh  (it  L2/ 2a )\  -  ( I  -  p) 


Ex  L  | 

2 p(\  P >  + 

yro 

cosh  (7 t  L:/2a) 

-  2/M  1  -  p) 

This  calculation  is  seen  at  line  6270  in  EETREN  It  reduces  to 


(26) 


.k>,„  =  rrr/MipV,,  Hid  -  s)/L  (27) 

at  line  6888  for  the  gradual  channel  case,  liquations  (25)  and  (27)  have  included  in  them  the  tempera¬ 
ture  correction  of  mobility;  also  Eq  (25)  has  in  it,  by  way  of  Eq.  (26),  the  correction  for  nonzero  chan¬ 
nel  thickness  in  Region  II  (picked  up  at  line  6325  where  the  hyperbolic  cosine  argument  is  adjusted). 


Resistance 


Output  Resistance 

In  addition  to  the  intrinsic  output  resistance  of  the  MESFET,  the  resistance  associated  with  the 
leakage  current  and  the  substrate  injection  current  must  be  considered.  It  is  probably  these  additional 
terms,  if  they  tire  disregarded  analytically,  which  make  the  measured  output  resistance  much  lower  than 
the  calculated  value  In  the  computer  program  EETREN  they  have  been  included. 


The  intrinsic  transistor  output  resistance  is  developed  analytically  as 


= 

9  Ins 


(28) 


In  St.it/  el  til  110)  the  development  of  the  expression  for  r0  is  concisely  shown,  for  the  assumption  of 
/ern-thickness-channei  in  Region  If  The  authors  did  the  full  derivation,  using  the  actual  channel 
thickness  Some  additional  terms  appear  in  the  expression  for  /•„  because  the  depletion  thickness  in 
Region  I!  is  a  function  of  //)s.  The  expression  from  Ref.  10  is: 


(/(/Pii.V,,  II  ( 1  />)  2 a 


2r'm 

"l|n|  <;<//<„  V„/;s  K 


TT  /.  i 

cosh  - 

2a 


1  1 
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The  lull  expression,  using  (a  8)  for  the  depletion  region  thickness  in  Region  II,  where  8  is  the  value 
given  by  Tip  (24),  consists  of  two  parts.  The  first  part  is 


urn 


|(a  SI 


(30) 


This  part  is  merely  Tig  (29)  with  (a  -8)  substituted  for  a.  The  second  part  is 


_ 2_ 

rr<//u  (A/)  W 


TT  L  2  IT  L  2  TT  L  2 

“sh  TTT^sT  '  smh 


(31) 


The  total  output  resistance  is  the  sum  of  r()|  and  r0  .  It  is  not  obvious  by  inspection,  but  after 
numerous  computer  runs  to  evaluate  the  significance  of  the  additional  part,  it  was  found  that 
almost  all  of  the  drain  resistance  correction  comes  from  the  substitution  of  (a  -  8)  for  a  in  Eq.  (29). 
Therefore,  in  FETREN,  Eq.  (30)  is  used  for  the  intrinsic  output  resistance  when  the  bias  allows  a 
Region  IT,  it  appears  at  line  6390  in  the  program.  When  there  is  no  Region  11,  then  L2  is  zero,  d 
replaces  p,  and  Li  is  L.  so  that  the  intrinsic  output  resistance  is  given  by  the  modified  second  term  of 
Eq.  (29): 


1  atinnNn  MM  1  -  </) 

which  is  the  same  as  that  developed  in  Van  de  Ziel  [71  and  Shockley  |3]  for  the  "gradual  channel"  case. 
The  calculation  of  Eq.  (32)  appears  in  the  program  at  line  6890 


Shunting  the  intrinsic  output  resistance  are  two  other  resistors.  The  first  is  to  account  for  simple 
surface  leakage,  and  substrate  leakage  due  to  overall  low  resistivity  and  lo  the  doping  tail,  as  was 
described  in  the  section  entitled  The  I-V  Curves  This  term  is  added  in  parallel  with  the  intrinsic  out¬ 
put  resistances  at  lines  6394  and  6896. 


The  second  shunting  resistor  is  related  to  the  substrate  injection  current  around  the  Gunn 
domain.  Assuming  that  most  of  the  drain  voltage  is  dropped  across  the  Gunn  domain  as  described  by 
Eastman  and  Shur  [9|.  then 


dlsu,  _  a  (Eq  23)  =  JT5 

6  *  />.s  a  *  ns  *  ns 


(33) 


The  reciprocal  of  this  is  the  resistance: 


(34) 

The  value  of  rmfl  is  calculated  in  FETREN  at  lines  6396  and  6892:  it  is  paralleled  with  the  intrinsic  out¬ 
put  resistance  at  lines  6398  and  6894. 

Parasitic  R csista nee 

The  parasitic  resistors  are  those  external  to  the  region  under  the  gate,  as  shown  in  Fig.  4.  The 
series  drain  resistor  and  source  resistor.  Rn  and  Rs  in  Fig.  4,  and  the  method  of  calculating  them  were 
described  in  the  first  paragraph  of  (he  section  entitled  The  Operating  Poini;  they  consist  of  bulk  sem¬ 
iconductor  resistance  and  ohmic  metal  contact  resistance. 

The  gate  resistor  (to  be  involved  in  noise  calculations  in  a  later  section)  is  a  bulk  resistance  of  the 
gate  metal  Fukui  (5 1  gives  a  suitable  method  of  calculating  its  value.  Ii  is  calculated  in  FETREN  at 
line  1370  and  used  at  line  9290  and  9300  in  the  noise  calculations. 


,=  2 


'St  H 
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Obviously.  if  ihe  values  of  these  parasitic  resistors  are  known  from  laboratory  measurements,  the 
measured  values  may  be  substituted  by  the  program  user  at  lines  1367,  1369,  and  1371 

I  he  tiu nn  domain  resistance  will  be  discussed  later,  in  the  section  entitled  Gunn  Domain. 

Input  Re\i\tance 

There  is  still  uncertainty  in  the  modeling  community  about  how  to  calculate  Rt\  of  Fig.  1.  The 
qualitative  form  of  R/N  versus  drain  voltage  for  a  particular  transistor  is  shown  in  Ref.  1 1  and  repeated 
here  among  the  measured  data  curves  to  be  discussed  later.  The  qualitative  form  of  Qv  versus  drain 
voltage  for  the  same  transistor  is  also  given  in  Ref.  i  I  and  is  repeated  here  among  the  measure  data 
curves.  It  is  assumed  that  for  drain  voltages  between  about  I  and  4  V  the  high  field  at  the  drain  end  of 
the  gate  is  in  the  transition  range  between  presence  and  absence  of  charge  accumulation  or  between 
presence  and  absence  of  a  dipole  layer  or  a  Gunn  domain.  This  range  is  difficult  to  model.  For  drain 
voltages  higher  than  about  3  V,  however,  an  inverse  proportionality  between  these  two  elements  is 
apparent.  Their  product  is  in  the  vicinity  of  the  transit  time.  So  in  FETREN  the  input  resistance.  RM , 
is  taken  as: 

*,v  =  CONST*  rw  •  (35) 

where  r  is  the  transit  time  to  be  discussed  in  the  section  entitled  Transit  Time.  The  resistance  is  calcu¬ 
lated  and  called  RINPHS  at  lines  6495  (if  there  is  a  Region  II)  and  7110  (if  there  is  no  Region  II)  for 
the  case  when  the  capacitance  calculation  disregards  the  Gunn  domain;  it  is  calculated  and  called 
RINSH  at  line  8662  if  a  Gunn  domain  is  considered.  In  the  interest  of  having  an  input  resistance  cal¬ 
culation  which  resembles  measured  data  for  low  drain  voltage,  the  transit  time  for  a  drain  voltage  of 
i'POl 2  is  used  when  the  drain  voltage  is  any  value  less  than  W2. 

Capacitance 


Input  Capacitance  or  Gate- ro-Source  Capacitance 


The  input  capacitance,  Qv.  has  two  different  names  (and  derivations)  in  FETREN;  it  is  CSG  at 
line  6491  and  7100  (no  Gunn  domain)  and  CGS  at  line  8631  (with  Gunn  domain).  In  both  cases  the 
capacitance  is  developed  as: 


However,  the  assumptions  made  to  determine  the  total  charge,  Q,  are  different. 


The  calculation  for  CSG,  assuming  no  Gunn  domain,  is  from  Pucel  el  al.  [4]  and  considers  only 
the  charge  in  Region  I  and  II  directly  under  the  Gate,  plus  a  fringing  term  which  is  independent  of  bias. 


The  calculation  for  CGS  is  from  Shur  112]  and  considers  the  charge  in:  (a)  the  depleted  area 
under  the  gate.  0V;  < b)  the  extended  depletion  area  beyond  the  gate  toward  the  drain.  0, ,  when  the 
high-field  domain  is  present;  and  (c)  the  depiction  area  in  the  high-field  domain  itself,  Q,hm,  when  it  is 
present.  Figure  6  shows  these  areas  conceptually.  The  theory  of  Shur  (12)  also  provides  a  valid  deriva¬ 
tion  of  the  input  capacitance  when  the  drain-to-source  voltage  approaches  zero;  so,  for  small  drain  vol¬ 
tage: 


/.  ft  4  ></  S  h 

tjt  7„;  t(; 


(37) 
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Fig  6  —  Conceptual  charge  arrangement  for  capacitance 
calculations  when  Gunn  domain  is  present 


(38) 


The  constant  PCDG  (set  at  line  8000  in  FETREN)  has  been  selected  for  best  agreement  with  available 
measured  data.  In  the  computer  program,  an  interpolation  routine  has  been  provided  to  calculate  a 
value  for  Q,  when  the  drain  voltage  falls  between  zero  and  half  the  pinch-off  voltage.  Experimental 
data  indicate  that  in  this  range  neither  Eq.  (37)  nor  (38)  predicts  the  capacitance  value  accurately.  The 
same  experimental  data  have  been  used  to  establish  the  interpolation  procedure,  performed  in  FETREN 
between  line  numbers  8710  and  8780. 


Feedback  Capacitance  or  Gate-to-Drain  Capacitance 

The  theory  of  Pucel  et  al.  (41  makes  no  provision  for  a  bias  dependent  feedback  capacitance ,  CFB 
in  Fig.  1.  It  is  there  treated  as  a  fixed  parasitic  capacitance,  which  will  be  described  in  the  next  subsec¬ 
tion.  Since  measured  data  show  CfB  to  be  bias  dependent,  in  general,  the  theory  of  Shur  [12]  is  used 
to  calculate  values  and  provide  them  as  CDG  in  FETREN.  For  the  same  reasons  as  were  outlined  in 
the  preceding  subsection,  two  ranges  of  validity  of  the  derivation  exist.  The  value  of  capacitance  which 
is  valid  for  drain  voltage  approaching  zero  is  equal  to  the  corresponding  gate-to-source  capacitance 
because  of  symmetry  at  zero  drain  voltage. 

C  =  C 
<*0  **0’ 

given  by  Eq.  (37).  The  value  for  large  drain  voltage  is: 

1/  1/  1/2  M  )1/2 

V o.  V, f  Afn 

CdK  =  1.46  (PCDG )  e0er  W  — -  —  . 

*ds  "<r 

Again,  available  experimental  data  have  been  used  to  interpolate  between  Eqs 
voltages  between  zero  and  half  the  pinch-off  voltage. 

Parasitic  Capacitances 

As  suggested  in  Pucel  el  al.  [4],  a  method  of  calculating  the  fixed  "interelectrode"  capacitances  is 
available  in  Smythe  [13],  with  some  interpretation.  Figure  7  gives  the  dimensions  which  are  used  in 
the  FETREN  calculations.  The  capacitance  between  any  two  electrodes  is  given  by: 

C  =  W  (tr  +  l)e0  K(JK\k:)k^  ■ 


(40) 

.  (39)  and  (40)  for  drain 
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Eig  7  —  Dimensions  for  inierelectrode  capacitance  calculations 


where  K  (k )  is  the  complete  elliptic  integral  of  the  first  kind  in  terms  of  the  modulus  k.  The  moduli  A, 
derived  from  Smythe  [131  for  the  relevant  geometries,  are  as  follows: 

1/2 

(42) 


(43) 

and 


Drain-Gate,  k<t  = 


Source-Gate.  A-  = 


(DL  +  DGL  +  WG)(DGL ) 


(DL  +  DGL)  (DGL  +  WG) 
( SL  +  SGL  +  WG)(SGL ) 


( SL  +  SGL)  (SGL  +  WG) 


Source-Drain,  ksJ 


(SL  +  DL  4-  SDL) (SDL) 
~(SL  +  SDL  )(DL  +  SDL) 


'44) 


In  FETREN  a  numerical  approximation  for  Ihe  complete  elliptic  integral  of  the  first  kind.  K  (k).  is 
included  as  a  subfunction  CEIFK,  calculated  at  lines  9100  to  9136.  This  was  done  in  keeping  with  the 
objective  ol  using  closed  form  solutions  throughout.  The  numerical  approximation  used  was: 


K(k)  =  [,4n-E  .4,(1  -  k2)  +  .4,(1  -  k2)2 1 


+  In 


[B0  +  B,<  I  -  k2)  +  BM  -  k2)2]. 


(45) 


where 


.4  n  =  I  3862944 
.4,  =  0.1119723 
.4,  =  0.0725296 
Bo  -  0.5 

=  0,1213478 
B:  =  0.0288729, 


15 


NKIDHRT  AND  SC'OTT 


giving  an  error  of  less  than  3  x  10  '  in  the  value  of  A' (A  ),  which  is  orders  of  magnitude  better  accu¬ 
racy  than  needed.  The  values  for  these  fixed  capacitances  are  calculated  between  lines  1315  and  1351. 
The  printout  for  them  is  at  line  1393.  Perhaps  the  most  significant  of  these  terms  is  CFSD,  the  output 
capacitance  (O0  in  Fig.  I). 


Sidewall  Capacitance 


Neither  theory  discussed  so  far  for  capacitance  calculations  (Shur  [12]  or  Pucel  et  al.  [4] >  provides 
for  the  extension  of  the  depletion  region  from  under  the  gate  in  the  direction  of  the  source.  In  Cooke 
114],  for  depletion  depth  less  than  the  full  active  layer  thickness,  the  charge  depleted  from  the  quarter- 
circle  adjacent  to  the  gate  is: 

Qsw  =  <1  Ndtt  -j-  W ;  (46) 


and  for  an  abrupt  junction  from  Poisson’s  equation,  the  voltage  is: 

2e  (fir 


so  that  the  capacitance  is 


C, 


stc 


dQsw 

dV 


(47) 


(48) 


where  z  is  the  radius  of  the  sidewall  depletion  area  (see  Fig.  3).  Also,  from  Cooke  [14],  when  the 
gate-to-source  voltage  is  larger  than  the  pinch-off  voltage,  an  approximate  value  for  the  sidewall  capaci¬ 
tance  is  given  by: 

cs"\r>rn-**r*  sin-'^F-’  (49) 

This  capacitance  term  is  calculated  in  the  computer  program  FETREN  between  lines  6500  and  6520.  It 
is  printed  out  as  the  third  term  in  line  6545.  This  capacitance  appears  as  C£X  in  Fig.  1.  Note  also, 
however,  that  the  interelectrode  capacitance  between  gate  and  source  metal  may  be  significant  in  C£X. 


Pad  Capacitance 


It  is  also  possible  to  have  significant  bond  pad-to-source  or  bond  pad-to-ground  capacitance. 
Cooke  (151  discusses  of  the  pad-to-source  capacitance.  It  consists  of  capacitance  to  the  substrate  offset 
by  the  Schottky  barrier  capacitance  formed  on  top  of  the  buffer  or  substrate  material.  The  pad-to- 
ground  capacitance  is  predominantly  a  "parallel  plate"  capacitance  which  may  be  significant  when  very 
thin  substrate  are  used.  Those  terms  are  process  dependent  and  may  not  be  significant;  they  are  not 
included  in  FETREN,  but  the  user  should  take  notice  of  their  possible  existence,  estimate  their  magni¬ 
tude,  and  include  them  separately  if  they  are  significant. 

Gunn  Domain 


As  described  earlier,  the  presence  of  a  stationary  dipole  domain,  or  a  Gunn  domain,  near  the 
drain  end  of  the  gate  has  a  profound  effect  on  the  input  capacitance  C//v,  the  input  resistance  R;v,  and 
the  feedback  capacitance  CfB.  If  a  Gunn  domain  can  exist,  there  is  no  indication  in  the  FETREN 
printout  — this  is  probably  the  most  common  case.  The  element  values  to  use  are  those  calculated 
assuming  the  presence  of  the  Gunn  domain  — namely,  CGS,  RINSH,  and  CDG.  Two  additional  terms 
are  calculated  and  printed  out  for  the  "Gunn  domain"  case:  CDOM,  the  domain  capacitance;  and 
ROOM,  the  negative  resistance  associated  with  the  domain,  both  discussed  in  Shur  [12].  They  are  cal¬ 
culated  at  lines  8460  and  8490  in  the  computer  program,  and  printed  out  at  line  8640.  These  terms  are 
shown  in  Fig.  1  as  CDOM  and  ROOM;  they  are  not  by  themselves  dominant  terms  in  the  equivalent 
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circuit— the  domain  effects  on  the  other  three  terms,  mentioned  above,  appear  to  be  much  more 
significant. 

If  a  Gunn  domain  cannot  exist,  the  printout  "NGD"  (No  Gunn  Domain)  appears.  In  this  case, 
the  element  values  to  use  are  those  calculated  assuming  no  Gunn  domain  — namely,  CSG  for  C/v, 
RINPHS  for  RIS,  and  CFDG  for 


The  calculation  for  C/v  and  ChB,  from  Shur  [12],  for  drain  voltage  approaching  zero  (CGSO  and 
CDGO,  respectively)  is  done  at  lines  8120  to  8140  and  printed  out  at  line  8180. 


The  criteria  discussed  in  Yamaguchi  et  al.  [16]  are  used  in  FETREN  to  determine  whether  a 
Gunn  domain  can  or  cannot  exist.  The  analysis  was  done  for  planar  geometry,  like  that  used  for  MES- 
FETs.  Poisson’s  equation  (Eq.  (50)),  the  current  continuity  equation  (Eq.  (51)),  and  the  current  den¬ 
sity  expression  (Eq.  (52))  were  solved  numerically  by  using  a  finite  difference  scheme: 

V20  -  (Nd  -  n),  (50) 


dti 

dt 


V  J, 


(51) 


J  =  nv  +  DV  n, 


(52) 


where  n  is  volume  electron  density,  t  is  time,  J  is  current  density,  v  is  drift  velocity,  and  D  is  diffusion 
coefficient.  The  result  relative  to  domain  formation  was  pleasantly  simple.  For  a  given  gate-to-source 
voltage,  the  depletion  layer  width  W,h ,  at  which  the  maximum  field  strength  just  equals  the  threshold 
field  Elh ,  for  bulk  negative  differential  mobility,  is  given  by: 


2e 


Wb,-  vgs  .  +  (£,/,  •  L)] 


1/2 


(53) 


Therefore,  if  the  total  active  layer  thickness  a  is  less  than  the  threshold  width  Wlln  a  Gunn  domain  can¬ 
not  form. 


The  calculation  of  W,h  is  at  line  number  8100  in  FETREN,  and  the  test  of  W,h  relative  to  the 
active  layer  thickness  is  at  line  8105.  If  a  Gunn  domain  can  exist,  the  program  proceeds  to  calculate 
the  relevant  parameter,  if  a  Gunn  domain  cannot  exist,  the  printout  "NGD"  appears  but  the  program 
proceeds,  anyway,  into  the  calculations  related  mostly  to  the  domain  case.  This  is  done  to  make  avail¬ 
able  to  the  program  user  the  value  for  CGSO  and  CDGO,  calculated  at  line  8137  and  printed  out  at  line 
8180;  and  the  value  of  transconductance  from  Shur  [12]  for  small  drain  voltage,  for  general  informa¬ 
tion.  This  value  of  transconductance  is  calculated  in  the  computer  program  as  BGM  at  line  8220,  and 
printed  out  at  line  8330. 

Incidentally,  if  the  active  layer  thickness  is  more  than  twice  Wlh,  the  Gunn  domain  can  even  pro¬ 
pagate.  In  summary  the  following  classification  can  be  made: 

a  <  tV„,  No  Gunn  domain, 

W,h  <  a  <  2  W,h  Stationary  Gunn  domain, 
a  >  2  W,h  Propagating  Gunn  domain. 

Active  layer  thickness  sufficient  for  domain  propagation  is  normally  not  seen  in  microwave  MESFETs; 
there  is  no  provision  in  FETREN  for  such  cases. 
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Transit  Time 


In  the  velocity  saturated  region  (Region  III  of  a  Ga  As  MfcSFET,  the  transit  time  is 

Z., 

T,  =  -  . 

v>„, 

In  Region  I, 

T\  =  - . 

V  a  vi' ra  tli’ 


(54) 


(55) 


where  an  average  velocity  can  be  determined  as 


^ averanv  ^  O^oirrani'  Mo 


Gn, 

it 


(56) 


and  where  I'/)S  =  TVo  (/>2  _  s2)  from  Ed-  (17),  and  L|  is  obtained  from  the  iterative  process  described 
earlier  and  has  the  value  indicated  by  Eq.  (12). 


If  there  is  no  region  II,  then  d  is  substituted  for  p  in  the  expression  for  rt,  and  there  is  no  t2. 
The  substitution  of  d  for  p  is  here  analogous  to  the  similar  substitution  in  the  drain  current  expressions 
(5)  and  (21). 


In  the  computer  program  FETREN,  the  calculations  for  transit  time  appear  at  line  6493  if  there  is 
a  Region  II,  and  at  line  7105  if  there  is  no  Region  II. 

NOISE  CALCULATIONS 

The  method  of  calculation  of  the  noise  parameters  is  taken  directly  from  Pucel  et  al.  14).  In 
FETREN,  printout  has  been  arranged  (at  line  9340)  for:  (a)  noise  figure  in  decibels  if  the  driving 
impedance  is  Z.S,  given  at  line  number  9156  (BFZSD);  (b)  noise  figure  in  decibels  for  optimum  noise 
driving  impedance  (BFMIND),  (c)  the  optimum  noise  driving  resistance  (BRSO);  and  (d)  the  optimum 
noise  driving  reactance  (BXSO). 

The  noise  calculations  appear  as  a  separate  subroutine  (ANF)  between  lines  9150  and  9343.  The 
gate  metal  thickness,  HGM,  is  assigned  at  line  944,  outside  the  subroutine. 

COMPARISON  WITH  MEASUREMENTS 


Figures  8  through  13  show  measured  and  calculated  values  of  the  predominantly  bias  dependent 
elements  of  a  particular  GaAs  MESFET.  The  measurements  were  done  at  NRL,  and  the  experimental 
results  were  reported  and  discussed  by  Willing  et  al.  [11).  The  calculated  values  are  from  FETREN, 
using  only  the  material  and  dimensional  parameters  given  in  Table  2. 


The  overall  agreement  between  measured  and  calculated  results  is  not  perfect,  but  it  is  reasonably 
self-consistent  and  is  better  than  any  other  method  known  over  the  full  range  of  both  of  the  bias  vol¬ 
tages.  No  noise  figure  measurements  were  made  on  this  device,  but  the  calculated  values  of  the  noise 
parameters  are  included  in  the  sample  printout  in  Appendix  A. 

The  values  for  the  fixed  resistances  and  capacitances  are  given  in  the  sample  printout  in  Appendix 
A  They  also  are  in  reasonable  agreement  with  measured  values. 

Calculations  and  measurements  of  various  parameters  have  been  compared  by  using  several  other 
transistors  for  more  limited,  special  cases.  Similar  fair  agreement  between  the  two  is  typical. 
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Fig  8(a)  —  Measured  static  characteristics 


Fig  8(b)  —  Calculated  static  /,,,-)', h  characteristics 
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VDS  (VOLTS) 


Fig.  9(a)  —  Measured  transconduclance 


VDS  (VOLTS) 


( ig  9(bl  —  Calculated  transconductance 


f 
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I  ig  10(b)  —  Calculated  inpul  capacitance 
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1.0  2.0  3  0  4.0  5.0 

VDS  (VOLTS) 

Fig  11(a)  —  Measured  feedback  capacitance 


I  tg  I  Kb)  —  C  alculated  feedback  capacitance 


(OR)  ohms  R„  ohms 
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Fig  13(a)  —  Measured  oulput  resistance 


Fig  13(h)  —  C  alculated  output  resistance 
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Table  2  —  Experimental  GaAs  MESFET  —  Physical  Parameters 


Gate  Width,  W 

600  p  m 

Gate  Length,  L 

1.7  p  m 

Active  Layer  Thickness,  a 

0.3 15  yum 

Donor  Impurity  Concentration,  .%> 

7.5  x  1  O' Vent' 

Low  Field  Mobility,  p0 

4240  cm2/(V  • 

Assumed  Saturation  Field,  t\ 

3.2  kV/cm 

High  Field  Saturation  Velocity  in  Gunn  Domain 

0.8  x  I07  cm/s 

Gate  Metal  Thickness,  HGM 

0.2  p  m 

GaAs  Relative  Dielectric  Constant,  er 

12.9 

Source  Metal  Length,  SL 

75  pm 

Drain  Metal  Length,  DL 

25  pm 

Drain-to-Gate  Spacing,  DGL 

2.15pm 

Source-to-Gate  Spacing,  SGL 

2.15pm 

Drain  Shunt  Resistance,  DR  I 

650  ohms 

Thermal  Resistance,  RTH 

1 17°C/W 

Characteristic  Doping  Density, 

3  x  lO'Vcm3 

Schottky  Gate  Built-in  Voltage, 

0.8  V 

Diffusion  Coefficient,  D 

45  ctn!/s 

CRYOGENIC  CONSIDERATIONS 

With  the  material  parameters  included  in  the  computer  program  FETREN,  it  would  seem  a  rela¬ 
tively  simple  matter  to  calculate  the  full  equivalent  circuit  and  the  noise  performance  of  a  MESFET  at 
any  temperature.  It  should  only  be  necessary  to  know  how  the  material  parameters  vary  with  tempera¬ 
ture.  This  may  indeed  be  true— and  it  may  not.  The  variation  of  the  noise  sources  themselves  with 
temperature  is  not  well  established;  also  the  difficulty  and  inaccuracy  associated  with  cryogenic  measure¬ 
ments  on  MESFETs  complicates  correlation  with  calculations,  especially  if  tuning  must  be  done,  with 
the  tuners  themselves  at  cryogenic  temperatures.  Weinrcb  1 17)  reported  that  the  correlation  between 
both  room  temperature  and  cryogenic  (20  K)  values  of  measured  minimum  noise  figure  were  in  good 
agreement  with  calculated  values  by  Pucel  et  al.  [4|,  the  method  used  in  FETREN;  but  at  the  same 
time  there  were  extremely  large  discrepancies  in  various  other  noise  coefficients,  raisings  questions 
about  the  validity  of  the  noise  figure  calculation  as  well.  This  area  needs  further  research 

DISCUSSION 

The  user  of  FETREN  should  be  familiar  with  FORTRAN  programming  and  editing.  Several  of 
the  input  parameters  must  be  written  in  by  the  user  before  running  the  program  The  inputting  pro¬ 
cedures  can  be  reprogrammed  by  the  user  to  suit  his  own  purposes. 

The  cost  of  the  sample  run  in  Appendix  A  was  a  bit  less  than  $8.  It  provided  154  calculated 
points,  al  a  cost  of  about  a  nickel  a  point. 

At  NRL  the  output  information  from  FETREN  has  been  used  as  input  information  for  the  well- 
known  time  domain  analysis  program,  SPICE.  Large-signal,  time  domain  analysis  of  a  2  GHz  MESFET 
amplifier  was  carried  out.  The  fundamental  output  signal  and  its  first  three  harmonics  showed  good 
agreement  with  measured  data.  This  was  so  in  spite  of  the  fact  that  the  FET  model  in  SPICE  is  the 
Schichman-llodgcs  model  of  a  junction  FET,  simulated  with  calculated  MESFET  values  from  FETREN. 

It  is  possible,  therefore,  to  go  from  semiconductor  material  parameters  and  dimensions  to  large 
signal  (as  well  as  small  signal,  low  noise)  performance  with  no  experimental  intermediate  steps.  This  is 
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irue  for  GaAs  MESFETs  at  or  near  room  temperature.  It  would  appear  to  be  a  fairly  simple  matter  to 

modify  the  program  FETREN  to  be  suitable  for  other  semiconductor  materials;  it  is  probably  also  possi¬ 
ble  to  modify  it  for  cryogenic  calculations,  but  neither  of  these  has  yet  been  done. 
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Program  FETREN  Listing,  Sample 
Printout,  and  Notation  List 
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100  PROGRAM  FETRFN ' INPUT • OUTPUTS 
150  GO  TO  917 
200  200  PRINT  300 

300  300  FORMAT '♦ENTER  A. CND. CND 1 . W6« U. H6M*> 

400  RFAP  ♦, A.CNP.CNm .U6.U.HGM 
500  500  PRINT  GOO 

GOO  GOO  FORMAT •♦FNTFP  IG.VGO.DVG»> 

700  READ  *. IG. VGO.DVG 
800  800  PRINT  900 

900  900  FORMAT  <*€NTFR  ID. VDO.  DVD*'* 

910  910  PRINT  915 

915  915  FORMAT  OFNTFR  A.CND»FO*> 

91G  READ  9.A.CND.FO 

917  917  A=. 315 

918  CND=7. 5E16 

919  FO=. 424 

922  CNDl=CNIi/l  0.  •♦lG 
924  CND2=CND1 
92G  TNOM=300. 

928  FO=FO*<  <TNOM.'300.  7>> 

930  UG= 1 . 7 
940  U=600. 

942  U1=U*.001 
944  HGM=. 2 
94G  UG1=WG*IE-4 
950  IG=3 
960  VGO=0. 

970  DVG=-2. 

1000  ID=3 
1010  VDO=l. 

,1020  DVD=2. 

1100  QEL=1 .602E-19 
1200  EPS0=8 . 854F- 1 2 
*  1280  EPS=12. 9 

1285  SEIF-0 

*  1290  CF=1 .56 

1300  3L=75. 

1303  DL=25. 

1306  DGL=2. 15 
1309  SGL=2.15 

'  1312  SDL=SGL*UG+DGt 

1315  CAPCO=EPSO*  < 1 . +EPS>  1 . EG 

1321  AKPG=$ORT  (  <  <PL  +DGL  +WG>  ♦PGL'*  ✓  <  <  DL+PGl  >  ♦  *PGL  +UF,>  >  ) 

1324  AKPDG=SGRT  < 1 . -AKDG*A* DG> 

*  1327  CFDG=CAPCa*CEIFK<A*PPG>/CEIFK  WHO) 

1333  ARSG=SORT  <  >.  <Sl  >SGL+UG>  ♦SGL>x  <  . SL+SGL  '  ♦  <SGL+WG>  >  > 

1336  AKPSG=SORT <1 . -AKSG*A* SG> 

1339  CFSG=CAPCO*CEIFK  <AKPSG>  -'CEIFK  <AKSG> 

1345  AKSP=30RT  <  <  (SI  +PI  *SPl  >  *SDL>  /  <  >  3L  +S DL  .*  ♦  <  DL  t-SDl-*  >  > 

1348  AKPSD=SQRT<1.-A*<SD*A*SD> 

1351  CFSD»CAPCO*CEIFK  <AKPSD> ^CEIF* <AFSD> 

1359  RDEN=W 1 *A*CNP 1 . 82 

1360  RC0=2.  \/ (Ul*A**.  5*CNP2**.66> 

1362  PCHS=1 . 1*S6L-  RDEN 

1364  RCHP= 1 . 1 *PGl /RPFN 
1366  RS=RCO+RCHS 
1368  RD=RCO*PCHP 

1370  PG= a 7. ♦ (W 1 ✓ 4. > **2 . > / <  WG*HGM*W 1 » 

1372  PRINT  1380  . RG.RS.RD 

1380  1380  FORMAT  (♦  (RG=*.F4.1»*  ;  RS=*. F4. 1 . ♦?  PP=*. F4.  1 , 

1393  PRINT  1396.CFPG.CFSG.CFSD 

1396  1396  FORMAT  o  <  CFDG=*.F7. 4> *PF  !  CFSG=*«F7. 4. *PF  !  CFSD=*.F7. 4. *PF 
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1397  PRINT  1398 

1398  1398  FORMAT '  * 

1400  VPO=OEL*A*fi*CNP*l .  E-6  '  <2.  ♦EPtO^EPT'1 
1410  CDR1  = 1 . 

1415  CDR2-3. 

1420  RSHUNT-650. 

1430  DR1=RSMLINT*<600.  ,'U> 

1450  RTH-70000.  ■'U 
1500  ES=3. 9 
1550  PRINT  1560 

1560  1560  FORMAT  <*  VG  VD  CSU  CUP 

1570  PRINT  1580 

1580  1580  FORMAT  (*  ♦> 

1605  F=FO 
1610  AVELSAT=.  8E5 
1615  0NCR-3.E 15 
1620  AK=4. 

1650  1650  PI=3. 1415927 
2000  XI  =  . 1*ES-*J6/VP0 
2100  VBI*. 8 
2200  DO  1100  1*1,16 
2300  VG*VGO+a-n4DV6 
2400  T=l. 

2600  USN-CVBI-VGi-'VPO 
2700  IF (WSN. GE. 0. )60  TO  1300 
2800  PRINT  1200 

2900  1200  FORMAT (*VG  GT  VBI  LINE  2600*> 

3000  RETURN 

3100  1300  S=SGRT<WSN) 

3150  ZP-S 

3200  IFCS.LE. l.iGO  TO  1500 
3300  PRINT  1400 

3400  1400  FORMAT OS  GT  1.0  LINE  3100*> 

3500  RETURN 

3600  1500  BET-2.  ♦A-'PI 

3800  DO  1100  J-t.ID 

3900  VD-VDO+-  <  J- 1 5  ♦DVD 

3901  3901  S=ZP 

3902  DET-1.5 
3918  3918  PL=S 
3920  PH=. 99999 

3930  VELSAT=ES*F*1.E5 
3935  GO=A*OEL*F*CND 
3940  AIS=GO*VI»€S*.  1 
3950  1600  STEP=<PH-PL> '10. 

4000  IF  <STEP. GT . 1 . E-8>  GO  TO  1800 
4100  PRINT  1700 

4200  1700  FORMAT  < *STEP  LE  lE-8*> 

4300  RETURN 

4600  1800  DO  1900  8 =1,11 

4700  P-PL+  <¥— 1 > ♦STEP 

4800  FI  =P*P-S»S-  <P**3-S0  3>  »2 .  -  3 . 

4900  AL  l=GO*%l*VPO*F  l'-'AI  S'  <  1 .  -P> 

5000  DEL  1 =VPO*  <P*R-S*S> 

5030  XYZ=<WG-AL1)/BET 

5040  XYZ=XYZ'CDR1 

5060  XYA=EXP  <XYZ> -FXP  <-XYZ> 

5080  XYB-XYA/2. 

5090  XYB=XYB*CDR 1 

5100  DEL2=VPO*-BET*<XI  /UG>  ♦XYB 

5120  BEL2=DEL2*CDR2 

5200  DEL=VP-DEL 1-PEL2-AIS*< 1 . -P> ♦ <RS+RD> 

5300  IF 'ABS 'DEL  > .LT. 1 . E-4) GO  TO  2200 
5400  IFrpEL.GT. 0. )GO  TD  2100 
5500  1900  CONTINUE 
56.00  PRINT  2000 


GM 
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5700  £000  FOPMATt^NO  CNVERGENCE  IN  11  ITERATIONS*) 

5750  RETURN 

'■EiOO  2100  PL=P-STEP 

59 On  PH=P 

6000  GO  TO  1600 

6010  2200  S2=S0RT  •  <  VB I -VG+AI S*  <  1 . -P)  *RS  -  VPO) 

6020  IF 'ABS <$2-S> -LT. l.E-4>60  TO  6100 

6030  S=S2 

6040  GO  TO  3918 

6100  6100  I F  <Ai 1 .  GT , UG>  GO  TO  6620 

6150  AL2=WG-AL1 

6151  BL21=AL2-'AL1 
6200  CUP=AIS*<1 . -P) 

6203  CVP=VP-CUR* <PS+RP> 

6206  CSUBA=A*><  'CNP  1 . E 1  ?>  ♦*.  25)  ♦  «C  VI'**. 5  ■  *M  1 .  E6 
6209  CSU»=CSUBA*  < AVFL  SAT,' .  8E5>  ♦  <  ■  CNC  R./? .  E  15  ■  ♦♦ .  25  > 

6212  CUP=i;UR*CSUB 

6216  CUR=CUP+VP  DR1 

6220  IF <PET.GT.2.>GD  TO  6320 

6225  PELT =CUP*1 . E6-  ■ OEL*CNP*W*VELSAT) 

6230  PP=VP*CUR 

6235  CPR1  =  <  A-.  25*PEL  T)  --B 

624  0  F=FO*  >:  ( ■  TNO(1+RTH*PP>  TNOM  >  .  ?.<  ■ 

6245  B£T=BET-*1. 

6250  GO  TO  3918 

632  0  632 0  X=P I ♦ <  UG-At  1 •  2.  A 

6325  ;,=X/CBR1 

633 0  Y« < EXP <  X) +EXP  ( -X) ) -  2 . 

634  0  AN=  c  < 1 . -$> *Y> - <1 . -P) 

635 0  AP=  < 2 . *P*  < 1 . -P>  +X I ♦AC 1  WG> *Y 
6360  PENA=AP-2.*P*<1.-P> 

6370  GM=AIS*AN--DENA  VPO 
6390  PP«VPO*PEN«/AI S-'d. -P > 

6394  PR=  <^PR1  *PP)  /  <PR1  ♦PR) 

6396  PSU=2.  *CVP-'CSUB 
6398  PR= <PR*PSU) / <PP+PSU> 

6400  F1=P*P-$*S- <P**3-S**3> *2.3. 

6410  FG=AN'PENA 

6420  F2=  <  (P**3-S**3> *2.  ••  3.  >  -  <  <P**4-S**4>  '2.  > 

6430  FC1  A=  <2.  *P*P*  <  1 .  -P)  *  <  1 .  -P)  +F2>  >•  <1.-P> 

6440  FC1B=S*  < 1 . -S) 

645  0  FC 1  =2.  *AL  1  ♦  <FG*FC  1 A-FC 1 B)  -'F 1  /  A 
6460  FC2A*2. ♦ (UG-AL 1 ) *FG/A 
6470  FC2B=1 . -2. *P*FG 

6480  2= (EXP <X> -EXP  < -X) ) /  <EXP  <X>  +EXP  <-X> ) 

6485  FC2C=<2. *UG*P'A/XI /Y>+2 

6490  FC2=FC2A+FC2B*FC2C 

649 1  CSG=EPS*EPSO*U* (FC 1 ♦FC2+CF) ♦ 1 . E6 

6492  IF'CSG.LT. <CFSG/2. > * CSG=CFSG^ 2. 

6493  TAU=  <Ai.  1  #AL  1  /  <F*VPO*  <P*P-S*S>  >  )  +1  0-  ♦  <UG-AL  1  >  /VELSAT 

6494  BFC=FC 1 +FC2+CF 

6495  PINPHS=TAU/CSG 

6496  IF  CVD.  EO.  0. ) RINPHS=TAU/CFSG 

6498  6498  IF <SE IF. EQ. 1 0) GO  TO  6560 

6499  IF  <SEIF.E0.20)GO  TO  6560 

6500  6500  CVGP=VB1-VG+CUP*PS 
6505  IF  <CVGP. GT . VPO) GO  TO  6520 
6510  CSU=PI*€PS*EPS0*U€*l.Ee/'2. 

6515  GO  TO  6545 

6518  6520  CSWl*SOPT'VPO  CVGP> 

6519  CSW2=AS IN  <CSW1 > 

6520  CSU=EPS*EPSO*WG*l . E8 

6521  CSW=CS«*CSW2 

6545  6545  PRINT  6550* VG.Vp, CSW- CUP. 6N» PR. CSG 

6550  6550  FORMAT <F6. 3. 2X. F7. 4. 2X, F6. 3*  3X. F5. 4 >  31  - F6 . 4 . 1 X. F8. 1 • 2X. E 1 0.  3> 
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6560  6560  60  TO  800O 
6570  6570  PBCD*3. 

6580  PRINT  7130.TPU.R1NPHS 

6585  CPU.  PNF  <P.S.F1 . PENP. BL81 » F6. X I . ft. PL 8* PI »F •OEL • WG«  Y.  CGS • Rt  «RG»  Gft. 
6586+THON . BFC . P I NSH> 

6590  SEIF=0 
6600  1100  CONTINUE 
6610  PETUPN 

6680  6680  CURL=0. 

6681  PL8=0. 

6688  BL81=0. 

6685  CUPH=PIS 

663 0  8  STEP=  t CUPH-C  URL  >  1 0 . 

6635  IF  <PB$  <  I.TEP>  .  6T.  1 .  E-8>  60  TO  3 
6637  GO  TO  65 on 
6640  PPINT  7 

6645  7  FDPMPT .♦5TEP  LE  1E-8.N0.8*> 

6650  PE TURN 

6655  3  DO  4  L  =  1 • 1 1 

6660  CUR=CURL>  <L-1>  *STEP 

6665  SSG=  <VBI-VG+CUR4PC :•  VPO 

6670  DSO=  <  VB I -V6+VP-C  UP*PP>  .  VPO 

6675  F8D=PS0-SSO-8.  ♦<prp**l .  5-?:'C'**l .  5> 

6680  E 1 =CUR-GQ*W*VP0*F8P-  WG 

6685  IF  tRBS ■ El> . LT. 1 . E-4 ■ GO  TO  5 

6690  IF  <E 1 . GE . 0 . > GO  TO  6 

6700  4  CONTINUE 

6708  GO  TO  6500 

6705  6  CURL  =CUP-STFP 

6710  CUPH=CUP 

6715  GO  TD  8 

6780  5  S»SORT<SSP> 

6785  P=SOPT • PSO ■ 

6830  FlD=P*Ii-S4$- '8.  3. • ♦ ■ D**3-:**3 . 

68  35  CUP*GO*W4VPp«F 1  P  UG 
6840  CVP=VP-CUP*  <PS+PP ■ 

6845  CSUBP=PK ♦ i <  C ND  1 . E 1 7 1 ♦« . 85 • ♦ • C VV*+. 5 ■ ♦U  1.E6 
6950  CSUB=CSUBft*<ftVELSPT  .  8E5  -  ♦  •  •CNCR  3.  G  15  ■  85  • 

686  0  CUP=CUP+C SUB 

6868  CUP=CUP+VP  PP1 

6964  IF • PET . GT . 8. > GO  TO  6888 

6868  PP=VP*CUP 

6878  F=E04<  ‘  <  TNOM+PTH*PP  1  TNON  7>  ■ 

6876  PET=PET+ 1 . 

6880  GO  TO  3918 
6888  6888  GM=GO*u* ■ P- I >  I.IG 
6890  pp=UG  U  GO  -l.-P1 
6898  PSU=8. *C VP  CSUB 
68*4  PP=  ■  PP*PT l.l .  •  PP+-PS  U  ■ 

68*6  pp= . pp i ♦pp .  . pp i +pp , 

6900  FG=  1  1  1 .  -S  ’  -  '  1 .  -P  1  >  '/I 

6*50  F8=  *  ■  v**?- :**'■< #8.  3.  ■  - 1 .  ti»54-:^^4  •  •  8.  * 

7o no  Fr  tP= 1 8. ♦P*P5 • 1 . -P ■ ♦ ■ 1 . -p> +F8 •  -  • 1 . -p ■ 

7  05  0  Fi;  t=  'FG»FC  IP-  :  ♦  ■  1  .  - 1  •  •  »8 .  ♦WG  F  1  P  ft 
'06  0  F:FC=FCl*i:F 

7!  0  0  r  :g=EPS*EPSD*i»#  •  Fr  1+CF  •  M  .  E6 
*  1  i  *4"  I F  ’  r  G .  L  T  •  *  l.  F  .  G  c'  .  1  1 1.  .  G=r  F  G  c" . 

7in5  Tpii=wG#uG  •F#vpp#.p^p-:^: .  ■ 

7iio  pinph'=tpu  csg 

7115  IF •VD.EO. 0. ■PlNPHi*Tftu  CFSG 

71  30  7130  FOPNPT.^  T  ftij®*  •  F8 .  8  •  *F  :EC  J  8  1 UFG  -  ♦,  F*.  8  .  ♦OHM  ;  ♦  > 

7800  GO  TO  64*9 
80  0  0  ft 0  00  Fi;PG=.  ;:8 

9040  pe : =  PVEL : ft? 
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X  MS  li  V'=I.I6*AE  ♦.  1 
SiTE.iT  AGQ=0EL4CNP*F ♦A#W  US 
807m  f  IJPFCH=AG04VS 
ft  08  it  PIF=45. 

8081  AV6=VG- ■  C  UP*P<T  1 
8 095  ETH=  •  4  385 .  -^Er  . E»F ■ *1 . E-4 

8100  UTH=-*  OPT  <  2 .  ♦FPt  0*FP<  ♦  •  Vf!  I  -  AVA  +  •  F  TH4MG  >  •  •  l!' FI  *C  HP  '  •♦lOOi  . 

810*5  IF  •  A-WTH •  81  10.81  10.8180 
81  1  0  81  10  PPINT  81  15 
8115  8115  FORMAT i*  NGP4> 

8180  8120  AD=SOPT  '2. ♦EPSC*EPS ♦ ' VB1 -AV6 ■  •  ©E  L  •  C  N  P  ■  ■  ♦  1  0 0 0 , 

81 22  AOP=AO*SOPT •  •VB1-VGT  < VP I -hVG >  ■ 

8i?5  ©=©el  ♦Ao*u*wE*rnri 

8130  AC  UP=A6C*  > 1 . -AO  A)  *V 3 

8135  CGSO=EPSC*EPS»W*WG*l .  E6-  ‘£.*f\0> 

8137  CGS0=C6S0*A0  APR 
8140  CBG0=C6S0 

8150  AGN=SOPT  (GEL  ♦CNP^EPS^EPSO  *  •  VP  I  -AVG  •  •  •  *F*€S  1 .  E5 

8180  PM=  1  0  0  0 . A6M 
8170  PM2=PM4*a. 

8175  IF  (T.GT.2, >G0  TO  8185 
8180  PRINT  8190.CGSO 

8190  8190  FORMAT  <♦  C6SO=CHGO=4* F7. 3* ♦FOR  VP=0. ♦> 

8195  8195  T=3. 

8200  AVI)=VP-CUR4 (RS+RB> 

8210  IF  <  RVP- VS > 822 0 . 822 0 . 837 0 

8220  8220  BGM=SQRT  <OEL*CNP*ERSO*€PS  1  2.  ♦  ■  VB I-AVG1)  )  >  «f«AVS*W/  (W6*T  .  E  -E  ) 

8230  PDQM=0. 

8240  BCUR=AG04  c  t .  -AO  A  -  *AVP 
8245  IFTSEIF.NF. 0>GC  TO  8637 
8320  PRINT  8330* BGM 

8330  8330  FORMAT (♦  BGM=* * F5 . 1  * ♦MMHO  ‘FOR  VPINT  LT  VS > ♦ ' 

8360  GO  TO  8637 

8370  8370  VPOM=AVP-VS 

8330  C=AGO*  <1 . -AO/-AT  *VS 

8460  C DOM- 364.  ♦  ‘A- AO)  *M#S&RT  (EP$4€RSO*OEL*SGRT  ‘CNCR^CNP  ■  VPOM  > 

8470  FM=. 00782«SGRT  (OEL^SGRT cCNB*CNCP  > ♦VPOM  <EPSO*EPS > • 

8480  DDOM*2060.  ♦SORT  ( VPDM*ERS*EPSO  >  ©EC ♦SORT  (CNIi^CNCR '  >  » 

8490  RDOM=l . 46E7*SORT  (EPS*EPSO»VPDM»4?>  (PIF*W* • A-AO ■ ♦ > OEL^OEL^CNP^CNCF  > *♦. 

8500  RV=<VBI-AVG>' VPOM 

8600  RN=SORT  (  SORT  CCNB< NCR>  ) 

8610  CDG=PCDG4l.46*EPS4€PSO*il4S©RT(PV)*RN*l .E6 
8620  CGS2=CPG-CPGxPV 
8625  CGS2=1 .64CGS2 

8630  CGS=ABS<-CGSO-eCGS2> 

8631  CGS=1 . 3*C6S 
8637  8637  GO  TO  870 n 

8639  IF  (SEIF. ED. 1(0  GO  TO  8690 

8640  SE40  PRINT  865  0. C  GS  *  C  BG . C  POM  *  ROOM 

8650  8650  FORMAT (•  CGS « C PG « C DOM « P P0M= ♦ « F6 . 3  *  F6 . 3 « P 7 . 5 « F8 . 1 ' 

8652  IF>VD.GE.VPOT)GO  TO  8660 

8654  PINSH=TAIJP  CGS 

C-E-S'E.  60  TO  t'EEc 

3660  8660  RINSH=TAU, CGS 

8662  8662  P IN"H=P IN "A  I  .  s 

8665  IF ■ VD.EO. 0. )PINSH=TAU  CFSG 

8670  PP1NT  8680* RINSH 

8630  3680  FORMAT.*  P1n:h=*»F7. 2 ■ 

3690  8690  GO  TO  6570 
8700  8700  I F • J . NE . 1 ■  GO  TO  8760 
8702  IF ■ 'FIF.FG.20.GO  Tn  8760 
8705  IF  • SEIF.EO. Id'  GO  TO  8740 

8710  ‘ E I F  = 1 0 
8715  VTEMP=VP 
8730  VPOT  =■  VPO  2 
8732  VPFI=VPO  7. 
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87  5? 

VP=VPOT 

8  7  •  7 

GO  TO 

39  0 1 

8740 

874  0 

SLOP  1 =  ■  C 

DG- 

»:  DGO 

•  VP 

874  1 

"L0P3 

= • CGSP  ■ 

-c\ 

.  .  V 

DF1 

8748 

SLOP  3 

=  *  .'GS-  '  C 

6  SO 

•  ■  C  '  * 

•  VPO*3 . 5 

874  5 

Y I  NT = 

• CGSO -3. 

:•  -  • 

SLOP 

3*VDF  1  ■ 

874*. 

*FIF  = 

8  0 

87?  0 

VD=VTENP 

87?8 

T8UP= 

TRIJ 

87?? 

GO  TO 

3901 

8760 

878  0 

IF'VD.GE 

.  VPOT  .  G 

0  TO  8640 

876? 

CDG=S 

LOP  1  *VD+i 

C  DG 

0 

8768 

IF  «  VD 

.  GE. VDF1 

■  iiO 

TO 

8778 

3770  rGS=3L0P3*VB+CGS0 

8771  GO  TO  8840 

8773  8773  CGS =SL0P3*VD+Y INT 

8730  GO  TO  8840 

9000  END 

9100  FUNCTION  CEIFK'RKO 
9103  80=1.3863944 

9108  Rl  =  . 1119733 

9109  83= . 0735398 
9113  B0=. 5 

9115  Bl=. 1313478 
9118  B3=. 0388739 
9131  RKC  1  =  1.  -RKC*fHf  C 
9134  RKC3=RKC  1  *W<C  1 

9137  RKC3=R 0+81  ♦RKC 1  +R8^RK C8 

9130  RKC4= (B0+B1 ♦RFC 1 +B3^RK C3> ♦RL OG ■ 1 .  88 C 1 > 

9133  CEIFK=RKC3+RKC4 

9138  RETURN 

9139  END 

9150  SUBROUTINE  8NF *P, S,F1 . PENR. BL31  • FG.XI, 8. 8L3«PI «F«OEL > WG« Y» CGf . PS.RG- 
9 1 5 1 *GN « TNON «  BFC  *RI NSH> 

9153  CONPLEX  2S- B2C » BY 1  IN -  BY 1 1 D.  BY 1 1 » BY31 . BRN 1 . BGN3« BZCR- BZC I . BZC  C 

9154  CGS=CGS»1.E-13 
9158  ZS=  <50. »  0. > 

9158  BF=1 0. E9 

9159  BW=3.*PI*BF 

9180  BP0=  <P*P-S*S- <  <P**3-S**3>  *4.  ,3.  >  ♦  •  .rp**4-s*#4'»  >  -FI 

9185  BDEL= I . 3 

9170  BPD=  <3. ♦BDEL/F1> ♦<  <  1 .  -P>  ♦♦3> ♦ < S-P+RLOG  •  •  1 .  -S  >  •  •  1 .  -P>  '■  > 

9173  BFR=PENR  <l.-p> 

9175  BG8M=  <  1 .  -P  >  ♦  <  1 .  -P)  ♦BFR-'F  1  Y 

9180  BK  =  <  <  <P**4-$*-*4  '  ,  8.  >  -  <  <:P**3-S#*S>  3.  >  +  (P~r  >  *  <  ' '  **d>  -  •  ,  ♦,=' .  ■>  p  ] 

9181  BK=BK>BG8N»P 

9185  BKP=BK ♦  5 1 . -p>  ♦  BL  3 1 

9190  BP08=BKP*BKP*  5P»P-S4S> -BKP*(BKP+BG8N> ♦ .p**  3-S»*3 ■ *4.  3. 

9191  BROB=0. 5*<FKP#FKP+4. ♦BKP»BG8N+BGRN*BGRM>  * • P**4-S **4 ■ 

9193  BPOC= < BG8N*BGRM* s 3. > -0. 8* iBKP^BGRN+BGRN^BGRM' * • P*»5 
9193  BPO=  < BROR+BROB+BROO x <F 1 1 ♦F 1 > 

9195  BPD8=  <3. ♦FKP^BGRN-BGRN^FGRN  >  *<P»P-S*S> 

9 1 98  BPDB=BGRM*BG8H*  <P++3-S+*3>  *S .  -- 3. 

9197  BPDC=3.  ♦  <  <BKP-BGRfU  ♦♦p >  ♦  <P-$-*-RLOG  •  •  1  . -P>  ✓  1 1 .  -S>  >  > 

9198  BPD=BDEL»v < I . -P> ♦♦S* ♦<BRPR-BRPB-BRpr >  *F1#F1#F1> 

9300  BS0R=  <  <P»P-S»S>  -  <P**3-S»*3>  *4.  -  3.  ♦  ■  P**4-S*»4>  ■'£.  )  ♦BKP 

9301  BSOB=  <  5P**4-S**4>  -  <  <P**3-S**3>  *£.  3.  >  -  >  P**5-S**5>  +£.  ..-'5.  >  ♦BGRN 
9303  BSO= <BS08*BS0B> / <F1#F1>  , 

9305  BDS8=  <  BK P-BG8N^ ♦  < S-P48L0G  <  <  1 .  -<"*>•  •  1 .  -P>  >  >  +BGRH*  <P*P-S*S>  '3 . 

9308  BSD=3. ♦BDEL*<  <  1. -P> ♦♦3) *BDS8-  'F1*F1 > 

93 10  BP  1  =  •  1  .  -P>  ♦  <  BPCI+PPP>  <F  1  •bG8N*BG8N*FG> 

9315  BP3R=P  I  ♦RLS'R 

9316  BP3B=  < EXP <BP38>  > - <4. *€XP  >  PP38  3. >  > *3. +BP38 

9317  BP3C=<1.-P>*PI'2. 

9318  BP3D=  'SIN'  BP3C  >  ■'  BP3C  >  **3 

9319  BDIF=35. 
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BPIF0=1 .  374F-83*TN0N*F*1  0000.  'QEL 
BF3=16. ♦BBJF*BP2P*BP2B  < <PI*43>  ♦BDIFCD 
BP2=  <  1 .  -P'  ♦WG*BF3'-  < X I *BFP*BFP*F6*P  > 

BKPP=BFP~B6PN*P  „ 

BP<?=  <  <WG  •  CR*XI ♦•  t  . -P'1  >  • *♦?'♦•  •  2.*F1*BF  PP -  BFP  BF C  '  **c  •  ♦Fo^l-F  -■ 
pp  j  p  j  ,  1  1  .  -P'1  9  >  ♦  '  c  .  ♦  I.H?  R  ' '  I  BGRM  BF  *_  ♦♦£  1  ♦F G*  •  B;F  0  +  BP  I'  1 

PL£TUG=RL2<  PIG 
BPP=BPl-*-BP2 
BPP=BP 1 ♦BP? 

BCCft-SGPT <  BP  1 ♦BP  1  BPP  BPP > 

BC C B  =  S OP  T  <  BP 2 ♦ BP2 - BPP  BP P  5 

BCCC=  1  BSO+BSP)  •  "OPT  <  1  PPO+-BPD  '  ♦  1  BPO+BPP>  > 

BC  C  =  PC  CC *BC  CP+BCCP 

BYt  1N=«;MPLX  •  0.  >  BW*CGS  > 

by 1 1  d=c riPL x  •  i .  -  B !>!♦•: g:*pin:h  ■ 

BY  1 1=BY1 IN  BYt IP 
BV£ 1 =GN  BY  1  IP 

BGGN=  ■  BW*BW*C  6"  ♦•!  G~  •  GM*BPP 
BGPN=GN*BPP 

BGN 1  =BY  1 1  BY£  1  OPT  1  BG ON  • 
pr,N<?=i:  MPL  X  «:  0 .  •  •  BC  C  ♦  -  OF  T  •  BGGN  •  •  ’ 

BGN=  1 C  ftB'"  1  BGN1  -BGN21  •  ♦♦2+  '  1 .  -BCC  ♦BCC  "•  ♦BGoN 
BPN=PG+PS+<BGPN  1  >  CRBi  '  BYc'l  ■  •  *♦£  •  ♦  •  1  .  -BC  c  ♦Be  C  1  ♦  BGGN  Bc-N  ■ 
BPN=BRN-  •  PG+PS  1  ♦  '  PG+P'""'  ♦TNOH  300. 

PrCG=BGPN^CPN  >G ■ BY  1 1  1  •  'CRB"  *BY21  '  •♦♦2* 

BZCB=BCC*SGPT  ■:  BGGN*BGPN '  BY£  1 
'  B7CC  = >  0.  *  1 .  • ♦BZCB 
i  B3C  =PG+P-  +  '  BZCR+BZCC  '  BgN 
,  BNP  1  =  >  C  RB “  •  Z  ‘ +Brc ■  1 ♦♦£  , 

;  BNFc  = ' BPN+BGN*BNF 1 >  • PERL  ■  2"  ■  • 

.  BF7'" P=  t  0 . ♦PL 06 t  0 • 1 . +BNF2 • 
i  BF :Ol=PERL ■ BZC ■ 

;  BP'0=C0PT ■ BPCD1**c+ ■ BPN  BGN • • 

I  Bx:0=-PINPG'B70 
>  BF  G=BF'P+BPP-2.  ♦BCC  ♦SOP T  <  BPP*BPP  ' 

(  BFC=  1  BPP-BCO^SOPT  •  BPP^BPP'*  '  BF  G 
;  BKP=BPP^BPP4- 1 .-BCC ♦BCC •  BF G 
i  BFMIN= 1 . ♦£. ♦BGN* • BPCOl+BPSD • 
r.  BFMINP=1  0.  ♦PLDGl  0  1  BENIN' 
i  PPINT  R. BEZCP«  BFNINP*  BP-0»  BXSO 

i  i  FppMftT  <  ♦  NZ=*'F*'.  1  •  ♦DB*  NM=*.F"?.  1  •♦OB!  POP=  ♦  <  F  .  .  1  >  ♦OHN  < 


1 • *00 


9343+N*> 

93*10  PPINT  9360 
9360  9360  FOPNRT  > ♦  ♦* 
9370  PETUPN 
9380  ENP 
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•  PG=  1.9  5  PS  =  4.0!  PP=  4.0' 

•  CFDG=  . 0677PF  !  CFSG*  . 0692PF  5  CFSIi=  .  094 IFF  • 

VG  VP  CSW  CUP  GM  DP  C  Zh 

0.000  1.0000  .030  .0893  .0388  3.8  . 346E-0J 

CGSO=C DGO=  . 47SFDP  VD=0. 

BGM=  32.1MI1H0  'FOP  VDIMT  LT  VS' 

CGS.CriG.CDOM.PDOn=  .883  .  38?  .  AlF.fi?  359.  1 

P I NS  H=  17.81 

TGU=  88.  65PSEC  !  PINPHS=  654. 74DHWS 

NZ=  5.3PB*  NM=  3.9PB!  POP=  33.00 HU!  XDP=  59.60HH 

0.000  3.0000  .030  .1458  .0543  101.8  . 769E+00 

CGS*  CBG» CB0M*PBC1H=  .881  .070  .01444  587.8 

PIHSH=  7.72 

TRU=  9. 12PSEC  !  PI NPHS*  1 l . 860HHS 

MZ=  9.3BB5  NH=  7.6DB!  POP=  18.50HH*  XOP=  31.50HH 

0.000  5.0000  .030  .1452  .0521  205.8  . 776E+00 

rGS. CDG. CPOM ■  PPOH=  .  754  .  045  .  0094 1  8  051.7 

PINSH=  5.88 

TRU=  8.40P3EC  !  P I NPHS=  10.S30HMS 

NZ=  9.6DB*  NM=  7. 7DB!  POP=  16.8DHM!  XDP  =  29.10HH 

-2.000  1.0000  .030  .0828  .0349  10.0  . 348E-01 

CGSO=CPg6=  . 852FOP  VD=0. 

CGS  * CPG f C DOM • P DOM=  .141  .  193  .02069  45.5 

PINSH=  31.81 

TGlJ=  1 3. 33PSEC  !  PI NPHS =  385.  30 OHMS 

NZ=  8. 4PB*  NM=  1.1DB!  PDP=  97.1DHM!  XDP  =  1 37. 30HM 

-8.000  3.0000  .030  .0738  .0408  165.9  .5G5E+00 

CGS » C DG . C DOM  *  F DOM=  .867  .086  .  00576  802  0.  1 

P IMSH=  16.30 

TRIJ=  8.28PSEC  !  PIHPHS  =  14.640HMS 

NZ=  4.6DB*  lift*  3.0IiB!  POP*  42.20HM*  XOF'=  77.0OHM 

-8.  000  5.0000  .090  .  0778  .  0  3*6  85,8.0  .567E  +  00 

CGS  *  C  DG*  CDOM • PDOM=  .386  .083  .00409  5518.7 

PINSH=  10.85 

TRIJ=  7. 55PSEC  5  PI NPHS*  13. 380HMS 

|SZ=  5.8DB!  NM=  3.6DB!  PDP=  32.90HM*  XOF'  =  k  1 . 50HM 

-4.000  1.0000  .  0  90  .  0088  .  0818  lM'.:-;  ,4-lE  +  oo 

MGD 

CGiO=CDGO=  . 193FDP  VD=0. 

NGD 

MGD 

C GS . C DG . f DOM *  PPOM=  .  101  .1 60  . 0 0 173  1972. 0 

PIMTH=  88.93 

TGI  l=  16.08PSEC  !  P I NPHS  =  37.280HMS 

MZ=  3.5DB!  NM=  1.5DB!  P0P=  188.0HHM!  XOP=  184.40HH 

-4.000  9.0000  .030  .0141  .0247  328.0  . 426E+00 

MGD 

ny-  *  c  DG,  r PPM .  PPC1M=  .154  .  O^G  .  00079  1810*2.5 

PINSH=  39.62 

TGI  1=  11. 57PSEC  !  PI  NPHS  =  27.140HM-: 

M7=  3.1PP!  NM=  1.5DB!  POP=  104.6OHM!  :'DP=  I42.80HM 

-4.000  5.0000  .030  .0187  . 0258  386.4  . 429E+00 

MGD 

CGS*CDG*CDOM.PPOM=  .239  .074  .0005.  44.25.8 

PINSH=  22.03 

TGU=  1 0. 02PSEC  !  P I MPHS  =  23. 36 OHMS 

MZ=  3.2DB!  NM=  1.9BB!  PQP=  70.6OHM!  XDP=  97.60HM 
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Program  FETRF.N  Notation  List 

In  pul  Para  met  e  r  s 

A  Active  layer  thickness,  /x  m 

CM)  Doping  level  under  gate/cm' 

CM) I  Doping  level  between  gate  and  drain  and  between  gate  and  source 

in  units  of  1 016/ cm' 

WG  Gate  length. /am 

W  Gate  width. /am 

HGM  Thickness  of  gale  metal,  /am 

1G  Number  of  gate  voltages  to  be  run,  integer 

VGO  Initial  gate  voltage,  V 

DVG  Gate  voltage  increment,  V 

ID  Number  of  drain  voltages  to  be  run,  integer 

VDO  Initial  drain  voltage,  V 

DVD  Drain  voltage  increment,  V 

FO  Low-field  mobility,  m2/(V  •  s) 

CND2  Not  used  in  program 

TNOM  Nominal,  ambient  operating  temperature,  K. 

QEL  Electronic  charge,  C 

EPSO  Free  space  permittivity,  e(),  F/m 

EPS  Dielectric  constant,  relative,  of  GaAs.  dimensionless 

SL  Source  metal  length,  /xm* 

DL  Drain  metal  length,  /xm* 

DGL  Gate  metul  to  drain  metal  spacing,  /xni* 

SGL  Gate  metal  to  source  metal  spacing,  /xm* 

SDL  Source  metal  to  drain  metal  spacing,  /xm* 

VPO  Pinch-oIT  voltage  (calculated) 

RTH  Assumed  thermal  resistance,  °C/W 

F.S  Saturation  field,  kV/cm 

AVELSAT  Saturation  velocity  for  extremely  high  fields,  m/s 
(for  Gunn  domain  calculations) 

CNC'R  The  characteristic  doping  level  (for  Gunn  domain  calcu!ations)/cni3 

VBI  The  gate  metal-to-semiconduvtor  Schottky  barrier  built-in  voltage,  V 

ZS  Complex  driving  impedance,  ohms,  for  calculation  of 

noise  figure,  NZ 

BF  Frequency  for  noise  calculations,  Hz 

BDIF  High  field  diffusion  coefficient.  cm2/s 

Output  Parameters 

Gate  metal  resistance,  ohms 

Gate-to-source  bulk  semiconductor  and  source  contact 
resistance  combined,  ohms 
Drain-to-gatc  bulk  semiconductor  and  drain  contact 
resistance  combined,  ohms 

Metal  pad-lo-metal  pad  fixed  capacitances.  pF 

Applied  gate  voltage,  V 

*Sw  ITT.  ? 


RG 

RS 

RD 

CFDG) 
CFSG  } 
C’FSD  ) 
VG 
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VD 

CSW 

CUR 

GM 

DR 

CSC 

CGSO 

CDGO 

BGM 

CGS 

CDG 

CDOM 

RDOM 

RINSH 

TAU 

RINPHS 

NZ 

NM 

ROP 

XOP 

NGD 


Applied  drain  voltage,  V 
Sidewall  capacitance,  pl; 

Drain-to-source  current,  mA 
Transconductance,  mntltos 
Druin-to-source  resistance,  ohms 
Source-to-gate  capacitance,  pP,  see  text 
Gate-to-source  capacitance,  pP,  when  VO  =  0 
Drain-to-gate  capacitance,  pl:,  when  VI)  =  0 
A  transconductance  calculation  by  1 12) 

Gate-to-source  capacitance,  pP,  see  text 
Drain-to-gate  capacitance,  pF 
Domain  capacitance,  pF 
Domain  negative  resistance,  ohms 
Input  resistance,  ohms,  see  text 
Transit  time,  ps 

Input  resistance,  ohms,  by  )4],  using  calculated 
transit  lime,  TAU 

Noise  figure,  dB,  when  driving  impedance  is  ZS 
Minimum  noise  figure,  dB,  for  optimum  driving  impedance 
Optimum  driving  resistance,  ohms,  for  minimum  noise  figure 
Optimum  driving  resistance,  ohms,  for  minimum  noise  figure 
"No  Gunn  Domain"  by  the  criteria  of  1 16] 


M  IDI  RT  AND  SCOTT 


Appendix  B 

Program  FETREN  Flow  Diagram 
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Appendix  C 

CALCULATION  OF  TEMPERATURE  EFFECTS  ON  GALLIUM  ARSENIDE  MESFETS 


R.  E.  Neidert 

Naval  Research  Laboratory,  Washington,  D.C.  20375 

ABSTRACT 

A  simple  way  of  accounting  for  temperature  effects  in  analytical,  equivalent  circuit  models  for 
microwave  gallium  arsenide  MESFETs  is  shown.  The  approach  is  based  on  the  knowledge  of  the  ther¬ 
mal  resistance  of  the  device  in  actual  operation,  and  the  variation  of  the  device’s  low-field  mobility  with 
temperature.  This  paper  develops  working  numerical  estimates  for  both  of  these  characteristics,  gives  a 
method  of  incorporating  them  into  the  model  calculations,  and  compares  measured  data  with 
temperature-corrected,  calculated  results. 

I.  INTRODUCTION 

Comparisons  between  measured  data  and  calculated  results  for  MESFETs  (metal-semiconductor 
field  effect  transistors)  might  agree  better  if  temperature  effects  were  included  in  the  calculations.  The 
method  of  accounting  for  these  effects  is  simple,  once  the  actual  thermal  resistance  of  the  device  in  its 
operational  configuration  is  known  and  a  relationship  between  the  device  temperature  and  its  low-field 
mobility  is  known. 

Section  II  provides  a  theoretical  basis  for  the  determination  of  the  above  two  unknowns  and  for 
their  utilization.  In  Section  III,  some  measured  and  calculated  results  are  compared. 

II.  THEORETICAL  CONSIDERATIONS 

The  variation  of  low-field  mobility’  with  temperature  is  a  strong  thermal  effect  in  microwave  GaAs 
MESFET  operation  [Cl].  The  rate  of  this  variation  is  itself  temperature  dependent  because  different 
carrier  scattering  mechanisms  predominate  in  different  temperature  ranges  [C2] :  also,  the  degree  of 
temperature  sensitivity  appears  to  be  doping-level-dependent,  as  is  shown  later.  For  purposes  of  this 
paper,  and  for  most  applications,  300  K  will  be  used  as  the  reference  environmental  temperature.  The 
method  used  to  develop  the  mobility-versus-temperature  relationship  centered  at  300  K  is  of  course 
applicable  to  any  other  reference  ambient  temperature.  (See  Ref.  C2  for  guidance.) 

The  published,  measured  Hall-mobility  data  given  in  Table  Cl  has  been  reviewed  in  an  effort  to 
find  a  suitable  mobility-versus-temperature  relationship.  The  distinction  between  Hall  and  Drift  mobil¬ 
ity  may  be  neglected  to  a  good  approximation  as  explained  in  Ref.  C3;  moreover,  only  the  slope  is 
involved  here.  Table  Cl  gives  the  exponent.  A,  relating  mobility  and  absolute  temperature,  in  the 
vicinity  of  300  K,  for  several  impurity  levels  of  n-type  GaAs,  according  to  the  relationship 

Mo11  T'. 

In  the  proportionality,  m  a  is  low-field  mobility  and  T  is  absolute  temperature  in  degrees  Kelvin.  The 
sources  of  the  data.  Refs.  C2  through  C6,  are  also  given  in  Table  Cl  A  linear  least-squares  fit  to  the 
first  live  data  points,  merged  with  the  asymptote  given  by  the  sixth  data  value  (for  intrinsic  material). 
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Tabled  —  Measured  Temperature  Sensitivity  of 
Mobility  versus  impurity  Level  (al  300  K) 


Reference 

Impurity  Level/cm' 

A 

~  C3,C4 

2.2  x  |(F 

-0.74 

C3,C5 

3.7  x  |015 

-0.97 

C2 

1.3  x  1015 

- 1 .46 

C2 

2.7  x  10lb 

-1.01 

C2 

7.1  x  1015 

-1.06 

C6 

intrinsic 

-1.80 

yields  the  curve  of  Fig.  Cl.  These  data  and  the  resultant  curve  are  certainly  not  conclusive;  but  the 
indication  is  that  in  the  vicinity  of  the  10,7/cnv'  doping  level  of  microwave  GaAs  MESFETs  an 
exponent,  -I,  of  -0.6  can  be  used.  In  the  convincing  comparison  with  experimental  results  given  in 
Section  111,  for  which  the  transistor  doping  level  is  7.5  *  I0"7cm’,  a  value  of  -0.7  is  used  for  the 
exponent  (as  suggested  by  Fig.  Cl).  Also,  in  subsequent  calculations,  the  saturated  velocity  was 
allowed  to  vary  directly  as  the  mobility,  as  was  suggested  in  Ref.  Cl. 


Fig  Cl  —  Fsliniaicd  ploi  of  mobility  temperature  sensitivity 
versus  impurity  level  at  300  K 


The  other  parameter  which  must  be  known  in  order  to  calculate  the  thermal  effects  is  the  device 
thermal  resistance  in  actual  operation.  This  value  can  be  measured  for  a  specific  case;  however,  it  can 
also  be  estimated  with  reasonable  accuracy.  All  device  manufacturers  are  aware  that  good  heat  sinking, 
especially  for  power  I  FTs,  is  important;  so  they  all  tend  toward  the  state  of  the  art  A  survey  of 
several  manufacturers  was  conducted  (in  a  few  cases  proprietary  information  was  provided),  and  pub¬ 
lished  information  was  reviewed,  to  determine  typical  thermal  resistances  The  values  were  found  to  be 
in  the  vicinity  of 

Rn,  =  (70000/  If  )  °C/W. 
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where  Rm  is  thermal  resistance  in  °C  per  watt  and  is  the  gate  width  in  microns.  This  equation  is  an 
oversimplification,  but  it  provides  a  useful  first  order  approximation.  It  is  possible,  of  course,  to  devi¬ 
ate  greatly  from  this  value  but,  based  on  the  findings  of  the  above  survey,  such  wide  deviations  are  sel¬ 
dom  found  in  practice,  either  for  power  or  low-noise  devices. 

The  above  information  forms  the  basis  for  calculating  the  device  characteristics  which  are  affected 
by  operating  temperature.  By  taking  advantage  of  the  fact  that  the  temperature  sensitivity  of  the  low- 
field  mobility  is  such  a  predominant  factor  in  overall  MESFET  temperature  sensitivity,  the 
temperature-corrected  values  for  the  dc  and  sma.I-signa!  characteristics  can  be  obtained  wilh  a  simple 
correction  in  the  mobility  value.  The  correction  procedure  consists  of  the  following  steps,  noting  that 
external  voltages  must  be  used  so  that  all  the  parasitic  resistances  are  included  in  the  power  dissipation 
calculation:  (a)  calculate  the  drain  current  at  th;  relevant  bias  point,  (b)  calculate  the  total  dissipated 
dc  power,  equal  to  drain  current  times  drain  voltage;  (c)  calculate  the  operating  temperature  of  the  dev¬ 
ice  by  using  the  total  dissipated  dc  power  and  the  thermal  resistance;  (d)  calculate  the  corrected  operat¬ 
ing  low-field  mobility;  (e)  repeat  steps  (a)  through  (d)  and  use  the  new  mobility  value  each  time  — until 
the  mobility  converges  to  its  final  value;  and  (f)  recalculate  all  characteristics  by  using  the  final, 
corrected  mobility  value. 

In  the  following  section,  the  mobility  correction  is  performed  at  each  bias  point,  for  the  experi¬ 
mental  device  doping  level  of  7.5  x  1016/cm\  according  to: 

Mo<at  T)  =  mo  (at  300  K) ((300  +  Rw  •  PD)/ 300)-°7, 
where  Pn  is  the  total  power  being  dissipated  by  the  device. 


(II.  COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED  RESULTS 

The  pertinent  information  on  one  of  the  transistors  used  for  experimental  comparison  is  given  in 
Table  C2.  The  calculation  method  was  that  of  |C7),  excep!  that  alternate  equations  were  derived  and 
utilized  when  the  method  of  (C7|  yielded  a  computed,  constant-mobility  region  length  under  the  gate 
which  was  longer  than  the  gate  itself.  (See  Refs  C 7,  C'8.  and  C 9  for  some  discussion  of  this  case.)  For 
temperature  correction,  the  mobility  was  adjusted  according  to  the  method  discussed  in  Section  II. 
Typically,  only  one  iteration  was  required. 


Table  C2  —  Parameters  of  Experimental  Tran-  istor 


Active  layer  thickness 

0  3 1 5  m  m 

Doping  level 

7  5  x  10"7cm' 

Gate  length 

1 .7  m  m 

Gate  width 

600  m  m 

Source  parasitic  resistance 

2.8  ohms 

Drain  parasitic  resistance 

3.7  ohms 

Intrinsic  GaAs  relative  diel.  const. 

12.9 

Thermal  resistance 

I17°C/W 

Low  field  mobility  at  300  K 

4240  cm2/V-s 

Schottky  gate  built  in  voltage 

0.7  V 
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The  calculated  and  measured  I-V  curves  of  the  device  are  given  in  Fig.  C2.  Two  sets  of  calculated 
curves  are  shown— one  with  temperature  correction  and  one  without.  The  similarity  between  the  meas¬ 
ured  curves  and  the  computed  curves  with  temperature  correction  is  convincing  evidence  of  the 
existence  of  strong  temperature  effects  and  the  suitability  of  the  method  of  accounting  for  them.  As 
would  be  expected,  temperature  effects  diminish  as  the  gate  voltage  approaches  pinch  off. 


Fig  02  —  Experimental  transistor  l-V  curves 


Without  showing  additional  curves,  the  following  discussion  outlines  some  other  observations  on 
temperature  effects  which  were  made  The  calculations  (both  with  and  without  temperature  correction) 
were  done  according  to  Ref.  C7.  Note  that,  for  constant  applied  gate  voltage,  an  increase  in  the  applied 
drain  voltage  results  in  an  increase  in  operating  temperature  The  temperature-corrected,  calculated 
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transconductance  decreased  with  increasing  drain  voltage  for  gate  voltages  in  the  vicinity  of  0  V,  con¬ 
sistent  with  measured  results.  On  the  other  hand,  the  equivalent  transconductance  values  calculated 
without  temperature  correction  increased  with  increasing  drain  voltage,  contradictory  to  the  measured 
data.  The  calculated  values  of  input  capacitance  were  nearly  unaffected  (less  than  2%)  by  the  tempera¬ 
ture  corrections.  Finally,  the  temperature-corrected,  calculated  drain  or  output  resistance  was  always 
higher  than  the  value  calculated  without  temperature  correction.  Therefore,  with  temperature  correc¬ 
tion,  the  increase  in  drain  resistance  with  increasing  drain  voltage  was  somewhat  larger  than  that  calcu¬ 
lated  without  temperature  correction. 

Fair  quantitative  agreement  was  obtained  between  the  measured  and  the  calculated  values  of  tran¬ 
sconductance  and  input  capacitance,  but  the  drain  resistance  calculated  by  Ref.  C 7  was  in  all  cases  much 
higher  than  the  measured  values  (approximately  four  to  ten  times  higher).  Therefore,  the  actual  sensi¬ 
tivity  of  the  output  resistance  to  temperature  is  not  necessarily  predicted  by  this  method  of  calculation. 
Work  is  continuing  to  include  additional  terms  which  contribute  to  the  output  resistance,  and  to  investi¬ 
gate  their  temperature  dependence.  The  additional  terms  being  studied  are  the  Gunn  domain  negative 
resistance  when  present  [C 10-02],  the  resistance  representing  hot  electron  injection  into  the  otherwise 
semi-insulating  substrate  when  the  Gunn  domain  is  present  [03],  and  simple  substrate  or  surface  leak¬ 
age  when  it  is  present. 

IV.  CONCLUDING  REMARKS 

The  effects  of  temperature  on  RF  performance  have  not  been  explicitly  presented  in  this  paper.  It 
is  possible  and  straightforward  to  calculate  these  RF  effects,  in  any  special  case,  by  computing  them 
from  the  modified  equivalent  circuit  element  values. 

A  second  and  somewhat  different  transistor,  having  300  nm  gate  width  and  1.0  /im  gate  length, 
with  an  active  layer  thickness  of  0.24  m  m,  has  also  shown  good  agreement  between  measured  I-V 
curves  and  calculated  curves  when  temperature  correction  was  used  in  the  calculations.  The  negative 
slope  in  the  I-V  curves  at  large  values  of  bias  power  appeared  again,  both  in  the  measured  data  and  in 
the  temperature-corrected  calculations. 

The  static  negative  resistance  in  the  I-V  curves  of  GaAs  MESFETs  is  not  necessarily  caused  by 
temperature  effects  alone.  There  can  be  other  effects  involved,  Norton  and  Hayes  [C14l  indicate  this 
with  some  reservations.  It  is  interesting  to  see,  however,  at  least  in  the  cases  considered  here,  the 
much  improved  agreement  in  the  calculated  I-V  curves  by  accounting  only  for  temperature  effects. 
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